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Abstract. 

We present JHKs observations of 22 intermediate-mass stars in the Scorpius-Centaurus OB association, obtained with the 
NAOS/CONICA system at the ESO Very Large Telescope. This survey was performed to determine the status of (sub)stellar 
candidate companions of Sco OB2 member stars of spectral type A and late-B. The distinction between companions and 
background stars is made on the basis of a comparison to isochrones and additional statistical arguments. We are sensitive to 
companions with an angular separation of 0.1" - 11" (13 - 1430 AU) and the detection limit is K s = 17 mag. We detect 62 
stellar components of which 18 turn out to be physical companions, 1 1 candidate companions, and 33 background stars. Three 
of the 18 confirmed companions were previously undocumented as such. The companion masses are in the range 0.03 M Q < 
M < 1.19 Mo, corresponding to mass ratios 0.06 < q < 0.55. We include in our sample a subset of 9 targets with multi- 
color ADONIS observations from Kouwenhoven et al. (2005). In the ADONIS survey secondaries with Ks < 12 mag were 
classified as companions; those with K s > 12 mag as background stars. The multi-color analysis in this paper demonstrates 
that the simple K s = 12 mag criterion correctly classifies the secondaries in ~ 80% of the cases. We reanalyse the total sample 
(i.e. NAOS/CONICA and ADONIS) and conclude that of the 176 secondaries, 25 are physical companions, 55 are candidate 
companions, and 96 are background stars. Although we are sensitive (and complete) to brown dwarf companions as faint as 
K s = 14 mag in the semi-major axis range 130 - 520 AU, we detect only one, corresponding to a brown dwarf companion 
fraction of 0.5 ± 0.5% (M i 30 Mj). However, the number of brown dwarfs is consistent with an extrapolation of the (stellar) 
companion mass distribution into the brown dwarf regime. This indicates that the physical mechanism for the formation of 
brown dwarf companions around intermediate mass stars is similar to that of stellar companions, and that the embryo ejection 
mechanism does not need to be invoked in order to explain the small number of brown dwarf companions among intermediate 
mass stars in the Sco OB2 association. 

Key words, binaries: visual - binaries: general - stars: formation - stars: low mass, brown dwarfs - associations - individual: 
Sco OB2 



1. Introduction 

The predominance of star formation in binary or multiple sys- 
tems inside stellar clusters makes the binarity and multiplic- 
ity of newly born stars one of the most sensitive probes of the 
process of star and star cluster formation (see Blaauw 1991, 
and references therein). Ideally one would like to have detailed 
knowledge of the binary population at the time that the stars are 
being formed. However, this is difficult to achieve in practice 
and therefore we have embarked on a project to characterize 
the observationally better accessible "primordial binary popu- 
lation", which is defined as the population of binaries as es- 
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tablished just after the gas has been removed from the forming 
system, i.e., when the stars can no longer accrete gas from their 
surroundings (Kouwenhoven et al. 2005). We chose to focus 
our efforts on the accurate characterization of the binary pop- 
ulation in nearby OB associations. The youth and low stellar 
density of OB associations ensure that their binary population 
is very similar to the primordial binary population. We refer 
to Kouwenhoven et al. (2005) and Kouwenhoven (2006) for a 
more extensive discussion and motivation of this project. 

Our initial efforts are concentrated on the Sco OB2 associ- 
ation. Sco OB2 consists of the three subgroups Upper Scorpius 
(US), Upper Centaurus Lupus (UCL), and Lower Centaurus 
Crux (LCC). The properties of the subgroups are listed in 
Table 1 . Its stellar population is accurately known down to late 
A-stars thanks to the Hipparcos catalogue (de Zeeuw et al. 
1999), and extensive literature data is available on its binary 
population (Brown 2001). In addition Sco OB2 has recently 
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US 
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64 


44 
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0.46 


0.67 


0.61 


UCL 


140 


15-22 


132 


65 


19 
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0.40 


0.61 


0.52 


LCC 


118 


17-23 


112 


57 
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0.37 


0.56 


0.44 


all 






308 


166 


36 
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0.41 


0.61 


0.51 



Table 1. Multiplicity among Hipparcos members of the three sub- 
groups of Sco OB2. The columns show the subgroup name (Upper 
Scorpius; Upper Centaurus Lupus; Lower Centaurus Crux), the dis- 
tance (see de Zeeuw et al. 1999), the age (de Geus et al. (1989); 
Preibisch et al. (2002) for US; Mamajek et al. (2002) for UCL and 
LCC), the number of known single stars, binary stars, triple systems 
and N > 3 systems, and the binary statistics (see § 7), after inclusion 
of the new results presented in this paper. 



been the target of an adaptive optics survey of its Hipparcos B- 
star members (Shatsky & Tokovinin 2002). We have conducted 
our own adaptive optics survey of 199 A-type and late-B type 
stars in this association (Kouwenhoven et al. 2005) using the 
ADONIS instrument, which was mounted on the ESO 3.6 me- 
ter telescope at La Silla, Chile. We performed these observa- 
tions in the K s -band (and for a subset of the targets addition- 
ally in the J and H band). We detected 151 stellar components 
other than the target stars and used a simple brightness criterion 
to separate background stars 1 from physical companions. All 
components fainter than Ks = 12 mag were considered back- 
ground stars; all brighter components were identified as candi- 
date companion stars (see also Shatsky & Tokovinin 2002). Of 
the 74 candidate physical companions 33 were known already 
and 41 were new candidate companions. 

In examining the binary properties of our sample of A 
and late B-stars we noticed that at small angular separations 
(< 4 arcsec) no companions fainter than Ks ~ 12 mag are 
present, assuming that the sources fainter than K s ~ 14 mag 
are background stars (as we had no information on their colors). 
The absence of companions with Ks > 12 mag and p < 4" is 
clearly visible in Figure 3 of Kouwenhoven et al. (2005). This 
result implies that A and B stars do not have close compan- 
ions with masses less than about 0.08 M Q , unless the assumed 
background stars are physical companions. In the latter case the 
close faint sources would be brown dwarfs (which are known 
to be present in Sco OB2; see Martin et al. 2004) and a gap 
would exist in the companion mass distribution. In either case 
a peculiar feature would be present in the mass distribution of 
companions which has to be explained by the binary formation 
history. 

We decided to carry out follow-up multi-color observations 
in order (1) to determine the reliability of our Ks = 12 mag 
criterion to separate companions and background stars, (2) to 
investigate the potential gap or lower limit of the companion 
mass distribution, and (3) to search for additional close and/or 
faint companions. 

These follow-up near-infrared observations were con- 
ducted with NAOS/CONICA (NACO) on the ESO Very Large 
Telescope at Paranal, Chile. We obtained JHKs photometric 

1 When mentioning "background star", we refer to any stellar object 
that does not belong to the system, including foreground stars. 



observations of 22 A and late-B members in Sco OB2 and 
their secondaries 2 . In Section 2 we describe our NACO sam- 
ple, the observations, the data reduction procedures, and pho- 
tometric accuracy of the observations. In Section 3 we describe 
the detection limit and completeness limit of the ADONIS 
and NACO observations. In Section 4 we determine the sta- 
tus (companion or background star) of the secondaries with 
multi-color observations. We perform the analysis for the sec- 
ondaries around the 22 NACO targets, and for those around the 

9 targets in the ADONIS sample for which we have multi-color 
observations. In Section 4 we also analyze the background star 
statistics, and we evaluate the accuracy of the K s = 12 mag 
separation criterion. In Section 5 we derive for each compan- 
ion its mass and mass ratio. In Section 6 we discuss the lack 
of brown dwarf companions with separations between 1" and 
4" (130-520 AU) in our sample, and discuss whether or not 
the brown dwarf desert exists for A and late-B type members 
of Sco OB2. Finally, we present updated binary statistics of the 
Sco OB2 association in Section 7 and summarize our results in 
Section 8. 

2. Observations and data reduction 

2.1. Definition of the NACO sample 

The major goals of our NACO follow-up observations are to 
determine the validity of the K s = 12 mag criterion that 
we used to separate companions and background stars in our 
ADONIS sample (Kouwenhoven et al. 2005), to study the com- 
panion mass distribution near the stellar-substellar boundary, 
and to search for additional faint and/or close companions. 

Our NACO sample consists of 22 member stars (listed in 
Table 2) in the Sco OB2 association: 10 of spectral type B, 

10 of spectral type A, and one each of spectral type F and G. 
The targets are more or less equally distributed over the three 
subgroups of Sco OB2: 9 in US, 6 in UCL, and 7 in LCC. 
All 22 targets are known to have secondaries in the ADONIS 
survey. 

We included in our sample all seven target stars with faint 
(Ks > 14 mag) and close (< 4 arcsec) candidate background 
stars: HIP61265, HIP67260, HIP73937, HIP78968, HIP79098, 
HIP79410, and HIP81949. The other 15 targets all have candi- 
date companion stars, for which we will use the multi-color 
data to further study their nature. Priority was given to tar- 
get stars with multiple secondaries (candidate companions and 
candidate background stars) and targets close to the Galactic 
plane (HIP59502, HIP60851, HIP80142, and HIP81972) be- 
cause of the larger probability of finding background stars. 

There are 9 targets in the ADONIS dataset with (photomet- 
ric) multi-color observations. These targets are also listed in 
Table 2 and all have secondaries. Kouwenhoven et al. (2005) 
use only the Ks magnitude to determine the status of a sec- 
ondary, including the secondaries around the 9 targets with 
JHK S observations. Later, in Section 4, we will combine the 
data of the 22 NACO targets and the 9 ADONIS targets with 

2 We use the term "secondary" for any stellar component in the field 
near the target star. A secondary can be a companion star or a back- 
ground star. 
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0.43 
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A3IV/V 


UCL 
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123445 
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0.41 


6.37 


B9V 


UCL 


73937 
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0.46 


6.23 


ApSi 


UCL 


78968 


144586 


5.87 


0.53 


7.42 


B9V 


US 


79098 


144844 


7.28 


0.45 


5.69 


B9V 


US 


79410 


145554 


7.00 


0.52 


7.09 


B9V 


US 


79739 


146285 


6.79 


0.52 


7.08 


B8V 


US 


79771 
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6.86 
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7.10 


B9V 


US 
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147001 


5.82 


0.44 
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UCL 
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5.80 


B5 V 


US 
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UCL 
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ADONIS multi-color subset 
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95324 


7.93 
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6.48 


B8IV 
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US 
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0.49 
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B9V 


us 


78530 


143567 
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0.48 


6.87 


B9V 


us 


78809 


144175 


7.20 


0.51 


7.51 


B9V 


us 


78956 


144569 


5.55 


0.50 


7.57 


B9.5V 


us 


79124 


144925 


6.41 


0.53 


7.13 


A0V 


us 


79156 


144981 


6.21 


0.53 


7.61 


A0V 


us 


80238 


147432 


7.64 


0.68 


7.34 


Al III/IV 


us 



Table 2. We have obtained follow-up multi-color observations with 
NACO for 22 targets in the Sco OB2 association. We include in our 
analysis 9 targets with multi-color observations in the ADONIS sam- 
ple. All targets listed above are known to have secondaries in the 
ADONIS survey. The table lists for each star the parallax and er- 
ror (taken from de Bruijne 1999), the K s magnitude, and the spec- 
tral type of the primary star. The Ks magnitudes are those derived 
in this paper for the 22 stars observed with NACO, and are taken 
from Kouwenhoven et al. (2005) for the other nine stars. The last 
column shows the subgroup membership of each star (US = Upper 
Scorpius; UCL = Upper Centaurus Lupus; LCC = Lower Centaurus 
Crux), taken from de Zeeuw et al. (1999). 



multi-color observations, and determine the status of the sec- 
ondaries of these 31 targets using their JHKs magnitudes. In 
the remaining part of Section 2 we describe the NACO obser- 
vations, data reduction procedures, and photometric accuracy. 

2.2. NACO observations 

The observations were performed using the NAOS/CONICA 
system, consisting of the near-infrared camera CONICA 
(Lenzen et al. 1998) and the adaptive optics system NAOS 
(Rousset et al. 2000). NAOS/CONICA is installed at the 
Nasmyth B focus of UT4 at the ESO Very Large Telescope 



on Paranal, Chile. The observations were carried out in 
Service Mode on the nights of April 6, April 28-30, May 4- 
5, June 8, June 19, June 25, June 27-28, July 3, July 24, and 
September 10, 2004. Some representative images are shown in 
Figure 1. 

The targets were imaged using the S13 camera, which has a 
pixel scale of 13.27 mas/pixel, and a field of view of 14" x 14". 
The CONICA detector was an Alladin 2 array in the period 
April 6 to May 5, 2004. The detector was replaced by an 
Alladin 3 array in May 2004, which was used for the remain- 
ing observations. We used the readout mode Double_RdRstRd 
and the detector mode HighDynamic. For both detectors, the 
rms readout noise was 46.2 e~ and the gain was « 1 le~/ADU. 
The full-well capacity of the Alladin 2 array is 4300 ADU, with 
a linearity limit at about 50% of this value. For the Alladin 3 
array the full-well capacity is 15 000 ADU, with the linearity 
limit at about two-thirds of this value. 

Each observation block corresponding to a science tar- 
get includes six observations. The object is observed with the 
three broad band filters J (1.253 //m), H (1.643 yum), and K s 
(2.154 /j.m). Since our targets are bright, several of them will 
saturate the detector, even with the shortest detector integra- 
tion time. For this reason we also obtained measurements in 
J, H, and Ks with the short-wavelength neutral density filter 
(hereafter NDF). The NDF transmissivity is about 1.4% in the 
near-infrared. The observations with NDF allow us to study 
the primary star and to obtain an accurate point spread func- 
tion (PSF), while the observations without NDF allow us to 
analyze the faint companions in detail. In order to character- 
ize the attenuation of the NDF we observed the standard stars 
GSPC S273-E and GSPC S708-D. These are LCO/Palomar 
NICMOS Photometric Standards (Persson et al. 1998). The 
near-infrared magnitudes of these stars are consistent with 
spectral types G8V and G1V, respectively. By comparing the 
detected JHK S fluxes, with and without the NDF, we deter- 
mined the attenuation of the NDF in the three filters. 

Each observation consisted of six sequential exposures of 
the form OSSOOS , where O is the object, and S is a sky ob- 
servation. Each exposure was jittered using a jitter box of 4 arc- 
sec. We used a sky offset of 15 arcsec, and selected a position 
angle such that no object was in the sky field. 

Each exposure consists of 5 to 35 short observations of in- 
tegration times in the range 0.35 - 5.3 seconds, depending on 
the brightness of the source. For the target observations without 
NDF we chose the minimum integration time of 0.35 seconds. 
For all standard star observations and target observations with 
NDF we chose the integration time such that the image does 
not saturate for reasonable Strehl ratios. We optimized the in- 
tegration time to obtain the desired signal-to-noise ratio. The 
short integrations are combined by taking the median value. 

Visual wavefront sensing was performed directly on the tar- 
get stars, which minimized the effects of anisoplanatism. For 
a subset of the target stars, the observation block was carried 
out multiple times. The target stars were usually positioned in 
the center of the field. Occasionally we observed the target off- 
center to be able to image a companion at large angular sepa- 
ration. 
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Fig. 1. With our NACO survey we find three close companions (shown in the figure), which were not detected in our ADONIS survey. The 
panels (6.6" x 6.6") are centered on the primary stars. Left: The binary HIP63204 in K$, with a companion at angular separation 0.15" and a 
background star at angular separation 1.87". Middle: The binary HIP73937 in K s , with a close companion atp = 0.24" and a background star 
atp = 3.56". Right: HIP79771 in K s , with two companion stars atp = 0.44" and p = 3.67". Companions are indicated with white arrows and 
background stars with black arrows. Several artifacts are visible in the fields of HIP63204 and HIP79771, which can easily be recognized as 
such. The panels show a subset of the total field of view for each observation, which is 14" x 14". For these three targets we observe no stellar 
components other than those shown in the panels. 



The observations on the nights of 29 April, 28 June, and 
3 July, 2004, were obtained under bad weather or instrumental 
conditions and were removed from the dataset. These obser- 
vations were repeated under better conditions later on in the 
observing run. The observations on the nights of 6 April and 
27 June, 2004, were partially obtained under non-photometric 
conditions, and were calibrated using the targets themselves 
(see § 2.5). All other observations were performed under pho- 
tometric conditions. Most observations (85%) were obtained 
with a seeing between 0.5 and 1.5 arcsec. For a large fraction 
of the remaining observations the seeing was between 1 .5 and 
2.0 arcsec. The majority (65%) of the observations were ob- 
tained at an airmass of less than 1.2, and for 98% of the obser- 
vations the airmass was less than 1.6. 

2.3. Data reduction procedures 

The primary data reduction was performed with the ECLIPSE 
package (Devillard 1997). Calibration observations, including 
dark images, flat field images, and standard star images, were 
provided by ESO Paranal. Twilight flat fields were used to cre- 
ate a pixel sensitivity map. For several observations, no twi- 
light flats were available. In these cases we used the lamp 
flats. The dark-subtracted observations were flatfielded and 
sky-subtracted. Finally, the three jittered object observations 
were combined. 

2.4. Component detection 

The component detection is performed with the STARFINDER 
package (Diolaiti et al. 2000). The PSF of the target star is ex- 
tracted from the background subtracted image. The flux of the 
primary star is the total flux of the extracted PSF. A scaled- 
down version of the PSF is compared to the other signals in the 
field with a peak flux larger than 2.5 - 3 times the background 



noise 3 . The profile of these signals is then cross-correlated with 
the PSF. Only those signals with a profile very similar to that of 
the PSF star (i.e., a correlation coefficient larger than « 0.7) are 
considered as real detections. Finally, the angular separation, 
the position angle, and the flux of the detected component are 
derived. 

As discussed in § 2.2, we observe each target in JHK$ with 
NDF to obtain an accurate PSF template, and without NDF to 
do accurate photometry on the faint companions. The observa- 
tions without NDF are often saturated, which makes PSF ex- 
traction impossible. None of the observations with NDF are 
saturated. Since these observations are carried out close in time 
and close in airmass, we assume that the PSFs of the obser- 
vations with and without NDF are not significantly different. 
This is illustrated in Figure 2, where we plot the radial profiles 
of the extracted PSF of HIP78968. For the saturated images we 
use the PSF that was extracted from the corresponding non- 
saturated image for analysis of the secondaries. 

Artifacts are present in the image when observing a bright 
object. The location of the artifacts is well-defined (CONICA 
manual), and since the artifacts look significantly non-circular, 
they can be easily recognized as non-stellar components (see 
Figure 1). Additionally, observations obtained with the NDF 
show a faint artifact at ~ 2 arcsec to the North-East of the target 
star. Care was taken that the extracted PSF and the analysis of 
companions in that area were not affected by the presence of 
this artifact. 



3 Note that a peak flux of 2.5 - 3 times the background noise corre- 
sponds to a total flux with a much larger significance, since the flux of 
a faint companion is spread out over many pixels. All detected com- 
ponents in our survey have a signal-to-noise ratio larger than 12.8. 
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Fig. 2. The radial profile of the PSF for the target star HIP78968. 
The observations are obtained using the NACO system in the night 
of May 5, 2004 in J, H, and K$. The corresponding Strehl ratios 
for these observations are 6.5% in J, 15.6% in H, and 23.6% in 
K s , Observations are obtained with neutral density filter (NDF; solid 
curves) and without NDF (dotted curves). The profiles are normalized 
such that the peak flux for the images obtained with NDF is 100 in J, 
10 in H, and 1 in K s . The observations with NDF allow us to extract 
the PSF and measure the flux of the primary star. The images obtained 
without NDF are much deeper, but the primary is usually saturated. 
Assuming the PSF is similar, we use the non-saturated PSF to analyze 
the secondaries in the saturated image. 

2.5. Photometry 

Observations of standard stars from the Persson et al. (1998) 
catalog were provided by ESO. The standard stars are used to 
determine the magnitude zeropoints for each night and each fil- 
ter individually. No standard stars are available in the nights of 
28 April and 10 September, 2004. For 28 April we determine 
the zero point magnitude using 2MASS as a reference system, 
and the target stars HIP67260, HIP67919, HIP68532 as sub- 
stitute standard stars. For 10 September, we use the 2MASS 
data for HIP83542 to calibrate J and H. For the K s filter we 
use the HIP83542 measurement of Kouwenhoven et al. (2005) 
since its 2MASS Kg magnitude is inaccurate due to confu- 
sion with the diffraction spike of a nearby star. The H mag- 
nitude of HIP59502 was obtained under non-photometric con- 
ditions on 6 April, 2004, so instead we used the corresponding 
2MASS measurement of this star. All observations in the night 
of 27 June, 2004 were obtained under non-photometric con- 
ditions; the fluxes of HIP80474 and HIP81972 are therefore 
calibrated using 2MASS. We derive the calibrated magnitudes 
using the mean extinction coefficients for Paranal: kj = 0.11, 
k H = 0.06, and k Ks = 0.07. 

The attenuation m NDF (in magnitudes) of the NDF is deter- 
mined for J, H, and K s using the standard stars GSPC S273- 



E and GSPC S708-D and the non-saturated target stars. All 
observations are performed in pairs (with and without NDF), 
which allows us to determine the NDF attenuation. For each 
filter m NDF is calculated as the median difference in magni- 
tude: m NDF = (m* >witho ut ndf - »J*,with ndf)- We find m NDF>J = 
4.66 ± 0.03 mag, m NDFiH = 4.75 ± 0.02 mag, and m NDF>Ks 
■1.85 + 0.02 mag, respectively. We do not measure a significant 
difference between the values found for the early-type program 
stars and the late-type standard stars. All observations done 
with the NDF are corrected with the values mentioned above. 

2.6. Photometric precision and accuracy of the NACO 
observations 

We estimate the photometric uncertainty of the NACO obser- 
vations using simulations (see § 2.6.1) and the comparison 
with other datasets (§ 2.6.2-2.6.4). For primaries, the exter- 
nal lcr error in J, H, and K s is ~ 0.04 mag, corresponding 
to an error of ~ 0.06 mag in the colors. Typical lcr exter- 
nal errors in magnitude and color for the bright companions 
(8 £ Kg /mag ^13) are 0.08 mag and 0.11 mag, respectively. 
For the faintest sources (Ks ^13 mag) the errors are 0.12 mag 
in magnitude and 0.17 mag in color. In the following subsec- 
tions we will discuss the analysis of our photometric errors. 

2.6.1. Algorithm precision 

The instrumental magnitudes of all objects are obtained 
using STAPvFINDER. We investigate the precision of the 
STARFINDER algorithm using simulated observations. We 
create simulations of single and binary systems with varying 
primary flux (10 4 - 10 6 counts), flux ratio (AK S =0-10 mag), 
angular separation (0" - 13"), Strehl ratio (1% - 50%), and 
position angle. We estimate the flux error by comparing the in- 
put flux with the flux measured by STARFINDER for several 
realizations. 

Using the simulations we find that the precision of the 
STARFINDER fluxes is ~ 1% (~ 0.01 mag) for most primary 
stars in our sample, for all relevant Strehl ratios and as long as 
the PSF of the primary star is not significantly influenced by the 
presence of a companion. The 1 % error is due to the tendency 
of STARFINDER to over-estimate the background underneath 
bright objects (Diolaiti et al. 2000). For the fainter primaries in 
our sample (flux between 5 X 10 4 and 5 x 10 5 counts), the flux 
error is 1 - 3% (~ 0.01 - 0.03 mag). 

For secondaries outside the PSF-halo of the primary, the 
error is typically ~ 0.01 mag if the flux difference with the 
primary is less than 5 mag. Fainter secondaries have a larger 
flux error, ranging from 0.01 mag to 0.1 mag, depending on the 
brightness of primary and companion. 

If the secondary is in the halo of the primary, its flux error 
is somewhat larger. For example, for a companion at p = 2" 
which is less than 4 magnitudes fainter than the primary, the 
flux error is 4% (~ 0.04 mag) or smaller. Deblending the PSF 
of a primary and close secondary does not introduce a much 
larger error, as long as the magnitude difference is less than 
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~ 5 magnitudes. No close companions with a magnitude differ- 
ence larger than 5 mag are detected in our NACO observations. 

For several fields the observations without NDF are satu- 
rated. In order to analyze the faint companions in the field we 
use the PSF of the corresponding non-saturated observation ob- 
tained with the NDF (see § 2.2). These observations are per- 
formed close in airmass and time, so that their PSFs are simi- 
lar. We estimate the flux error by comparing PSFs correspond- 
ing to non-saturated images obtained with NDF and without 
NDF. These comparisons show that the resulting error ranges 
from 0.02 to 0.5 magnitudes, depending on the brightness of 
the secondary. We therefore minimize flux calculations using 
this method, and only use measurements obtained with the PSF 
of the non-saturated image when no other measurements are 
available. In the latter case, we place a remark in Table A. 1 . 

2.6.2. Comparison with 2MASS 

We compare the near-infrared measurements of the 22 targets 
in our NACO survey with the measurements in 2MASS (Cutri 
et al. 2003) to get an estimate of the external errors. We only se- 
lect those measurements in 2MASS that are not flagged. Since 
the resolution in our observations is higher than the « 4" res- 
olution of 2MASS, we combine the observed fluxes of the 
primaries and close companions before the comparison with 
2MASS. For the observations not calibrated with the 2MASS 
measurements, the rms difference between our measurements 
and those of 2MASS are 0.055 mag in J, 0.040 mag in H, and 
0.049 maginlj. 

2.6.3. Comparison with the ADONIS survey of 
Kouwenhoven et al. (2005) 

We detect all but two of the stellar components found by 
Kouwenhoven et al. (2005) around the 22 target stars in our 
NACO survey. We do not observe the faint companions of 
HIP80142 atp = 8.54" and HIP81949 at p = 9.70" because 
they are not within our NACO field of view. We find three 
bright companions at small angular separation of HIP63204 
(p = 0.15"), HIP73937 (p = 0.34"), and HIP79771 (p = 
0.44"). Since these objects are not found in the ADONIS obser- 
vations of Kouwenhoven et al. (2005), their fluxes and those of 
the corresponding primaries are summed for comparison with 
Kouwenhoven et al. (2005). 

The rms difference between the 22 primaries observed with 
NACO and those described in Kouwenhoven et al. (2005) is 
0.055 mag in K s . HIP69113 and HIP78968 additionally have 
multi-color observations in Kouwenhoven et al. (2005), which 
are in good agreement with the measurements presented in this 
paper. The J and H measurements of HIP80474 and HIP80799 
are flagged "non-photometric" in Kouwenhoven et al. (2005), 
and are not discussed here. 

Our dataset and that of Kouwenhoven et al. (2005) have 35 
stellar components other than the target stars in common. The 
rms difference between the K s magnitude of these objects in 
the two papers is 0.26 magnitudes. The differences are similar 



for the objects that have common J and H measurements in 
both papers. 

2.6.4. Comparison with Shatsky & Tokovinin (2002) 

Three targets in our NACO survey are also included in the 
binarity survey amongst B-stars in Sco OB2 by Shatsky 
& Tokovinin (2002): HIP79098, HIP80142, and HIP81972. 
Seven secondaries are detected both their survey and in our 
NACO survey (1 for HIP79098; 2 for HIP80142; 4 for 
HIP81972). Shatsky & Tokovinin (2002) classify these seven 
secondaries all as 'definitely optical' or 'likely optical'. They 
performed their observations in both coronographic and non- 
coronographic mode, and were therefore able to find five faint 
secondaries which do not appear in our NACO sample. 

The J and K s magnitudes of HIP79098 and HIP80142 
and their companions are in good agreement with our mea- 
surements. Our measurements of HIP81972 are in good 
agreement with those in 2MASS as well as the measure- 
ments in Kouwenhoven et al. (2005), but there is a discrep- 
ancy between our measurements of HIP81972 and those in 
Shatsky & Tokovinin (2002). The magnitude difference be- 
tween HIP81972 and its companions in our observations and 
in Shatsky & Tokovinin (2002) are similar. The observations 
of HIP81972 are flagged 'likely photometric' in Shatsky & 
Tokovinin (2002), but since they disagree with those in this pa- 
per and those in 2MASS, we assume they are non-photometric, 
and ignore them for the magnitude comparison. 

2.7. General properties of the NACO observations 

In the fields around the 22 targets we observed with NACO, we 
find 62 components other than the target stars. The properties 
of these targets and their secondaries are listed in Table A. 1 . 
The 22 primaries have 5.5 mag < J < 7.8 mag, 5.0 mag < H < 
7.7 mag, and 4.9 mag < K s < 7.7 mag. The brightest com- 
panions observed are ~ 7.5 mag in the three filters, while the 
faintest secondaries found have J = 16.6 mag, H — 17.3 mag, 
and K s = 17.3 mag. 

With NACO we detect components in the angular separa- 
tion range 0.15" < p < 11.8". The lower limit on p depends 
on the Strehl ratio and the magnitude difference between pri- 
mary and companion (see also § 3). The upper limit is de- 
termined by the size of the field-of-view. The median formal 
error in angular separation is 4 mas for bright components 
(8 mag ^ K s ^ 13 mag) and 10 mas for faint components 
(Ks ^13 mag). Position angles are measured from North to 
East. The median formal error in the position angle is 0°.007. 

We find 27 stellar components that are not detected by 
Kouwenhoven et al. (2005). Three close secondaries are found 
at small angular separation from HIP63204-2 (p = 0.15"; 
K s = 8.40 mag), HIP73937-1 (p = 0.34"; K s = 8.37 mag), 
and HIP79771-2 (p = 0.44"; K s = 11.42 mag). The for- 
mer two have been reported as candidate companions (Worley 
& Douglass 1997); the latter was previously undocumented. 
The other 25 new secondaries are all faint (K s ^ 12 mag). 
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Fig. 3. The detection limit and completeness limit for several targets 
in our ADONIS survey (Kouwenhoven et al. 2005) and NACO sur- 
vey (this paper). The target star, Strehl ratio, and instrument are in- 
dicated in each panel. The lower and upper curve show the detection 
limit and the completeness limit, respectively. The detection and com- 
pleteness limits shown above are representative for the ADONIS and 
NACO observations. The target stars have K s magnitudes of 7.33 mag 
(HIP81949), 7.42 mag (HIP78968), 7.45 mag (HIP80799), 6.52 mag 
(HIP58859), 6.26 mag (HIP76048), and 6.71 mag (HIP65178). The 
completeness limit is generally ~ 0.3 mag brighter than the detec- 
tion limit. At close angular separation a higher Strehl ratio results in a 
fainter detection limit. The dashed rectangle encloses the region with 
12 mag < K s < 14 mag and 1" < p < 4", which is relevant for our 
analysis of the substellar population in Sco OB2 (see § 6). In this re- 
gion the NACO observations are complete and ADONIS observations 
are more than 95% complete. 

Two of these 25 secondaries were also reported by Shatsky & 
Tokovinin (2002). 

3. The completeness and detection limit of the 
ADONIS and NACO surveys 

We cannot detect sources fainter than a certain magnitude be- 
cause of the background noise in the images. The faintest de- 
tectable magnitude additionally depends on the angular dis- 



tance to the primary star, the primary star magnitude, and the 
Strehl ratio. For a correct interpretation of the results of the sur- 
vey it is therefore important to characterize the limiting magni- 
tude of the observations (the detection limit) and the magnitude 
at which a star is likely detected (the completeness limit). 

We study the completeness limit and detection limit as 
a function of angular separation from the primary for six 
stars. These are HIP58859, HIP65178, and HIP76048 from 
our ADONIS survey and HIP80799, HIP78968, and HIP81949 
from our NACO survey. These stars are selected to cover the 
range in Strehl ratio of the observations, so that the complete- 
ness and detection limits are representative for the other targets 
in the ADONIS and NACO surveys. 

For each observation STARFINDER extracts the PSF from 
the image (see § 2.4). We simulate observations by artificially 
adding a scaled and shifted copy of the PSF to the observed im- 
age. We reduce the simulated image as if it were a real observa- 
tion. We repeat this procedure twenty times for simulated sec- 
ondaries with different angular separation and magnitude. We 
define the detection limit and completeness limit as the mag- 
nitude (as a function of angular separation) at which respec- 
tively 50% and 90% of the simulated secondaries are detected. 
The curves in Figure 3 show the completeness and detection 
limit for the six stars mentioned above. Due to our sampling 
the magnitude error of the completeness and detection limit is 
~ 0.15 mag. The figure clearly shows that a high Strehl ratio 
facilitates the detection of closer and fainter objects as com- 
pared to observations with lower Strehl ratio. For all stars the 
completeness limit is ~ 0.3 mag above the detection limit. 

In Section 6 we analyze the substellar companion popula- 
tion in Sco OB2 in the angular separation range 1" < p < A" 
and magnitude range 12 mag < Kg < 14 mag. The NACO 
observations of 22 targets are complete in this region. Only a 
few targets in the ADONIS sample are incomplete in this re- 
gion. Assuming a flat semi-major axis distribution, we estimate 
that about 5% of the faint companions at small angular separa- 
tion are undetected for the 177 targets that are only observed 
with ADONIS (see Figure 3). For the combined NACO and 
ADONIS sample this means that we are more than 95% com- 
plete in the region 12 mag < K s < 14 mag and 1" < p < 4". 

With NACO we are sensitive down to brown dwarfs and 
massive planets. To estimate the mass corresponding to the 
faintest detectable magnitude as a function of angular separa- 
tion, we use the models of Chabrier et al. (2000). We assume a 
distance of 130 pc, the mean distance of Sco OB2, and an age 
of 5 Myr for the US subgroup and 20 Myr for the UCL and 
LCC subgroups (cf. § 4. 1). At a distance of 130 pc, the bright- 
est brown dwarfs have an apparent magnitude of K s « 12 mag. 
With NACO we are able to detect brown dwarfs at an angular 
separation larger than « 0.3". The magnitude of the faintest 
detectable brown dwarf increases with increasing angular sep- 
aration between p = 0.3" - 2". The minimum detectable mass 
is a function of age due to the cooling of the brown dwarfs. 
For angular separations larger than « 2" the halo of the pri- 
mary PSF plays a minor role. For these angular separations we 
are sensitive (but not complete) down to K s « 16.5 mag with 
NACO, corresponding to planetary masses possibly as low as 
~ 5 Mj for US and ~ 10 Mj for UCL and LCC. However, 
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because of the large number of background stars and the un- 
certain location of the isochrone for young brown dwarfs and 
planets, we will not attempt to identify planetary companions 
in this paper (see § 4.2 for a further discussion). 

Several secondaries are not detected in one or two fil- 
ters. For the filter(s) in which the secondary is not observed 
we determine a lower limit on the magnitude of the miss- 
ing star using simulations. We perform the simulations as de- 
scribed in Kouwenhoven et al. (2005). The position of the sec- 
ondary is known from observations in other filters. This po- 
sition is assigned to the simulated companion, so that the de- 
tection limit as a function of angular separation is taken into 
account. Images are created with simulated secondaries of dif- 
ferent magnitude. We reduce the images as if they were real 
observations. The lower limit on the magnitude is then deter- 
mined by the faintest detectable simulated secondary. Two sec- 
ondaries (HIP73937-1 and H1P76071-1) have a lower limit in 
J because they are unresolved in the wings of the primary star 
PSF. The other secondaries with a lower magnitude limit for a 
filter have a flux below the background noise. 

4. Status of the stellar components 

In this section we determine the status (companion star or back- 
ground star) of the secondaries. We analyze the secondaries de- 
tected around the 22 targets observed with NACO, as well as 
the secondaries around the 9 targets with multi-color observa- 
tions in the ADONIS sample. For the 31 targets analyzed in 
this section we detect 72 secondaries in total. Sco OB2 mem- 
bers and their companions should be located near the isochrone 
in the color-magnitude diagram, while the background stars 
should show a much larger spread. We use this property in § 4.2 
to separate companions and background stars. As a consistency 
check we study in § 4.3 how our results compare to the ex- 
pected number of background stars in our observations. 

4.1. Color-magnitude diagram and isochrones 

For the 22 targets in the NACO sample and the 9 targets in 
the ADONIS multi-color subset we have magnitudes in three 
filters, as well as for most of their secondaries. Several of the 
faintest secondaries are undetected in one or two filters, as they 
are below the detection limit. Color-color and color-magnitude 
diagrams with the 3 1 targets and 72 secondaries are shown in 
Figure 4. Lower limits are indicated for objects that are not 
detected in all three filters. Several secondaries in our sample 
are either very red or very blue. These secondaries fall outside 
the plots in Figure 4, and are all background stars. 

For our analysis we adopt the isochrones described in 
Kouwenhoven et al. (2005), which consist of models from 
Chabrier et al. (2000) for 0.02 M < M < 1 M , Palla & 
Stahler (1999) for 1 M < M < 2 M Q , and Girardi et al. (2002) 
for M > 2 M Q . For members of the US subgroup we use the 
5 Myr isochrone, and for UCL and LCC members we use the 
20 Myr isochrone. 

Absolute magnitudes Mj, Mh, and Mk s are derived from 
the apparent magnitudes using for each star individually the 
parallax and interstellar extinction Ay from de Bruijne (1999); 



see Kouwenhoven et al. (2005) for details. The error on the 
parallax is 5 - 10% for all targets, and can therefore be used to 
derive reliable distances to individual stars (e.g., Brown et al. 
1997). The median fractional error on the distance is 6.5%, 
which introduces an additional error of 0.15 mag on the ab- 
solute magnitudes. Combining this error with the error in the 
apparent magnitude (§ 2.6) we obtain lcr error estimates for the 
absolute magnitudes of 0.16 mag for the primaries, 0.17 mag 
for the bright companions, and 0.19 mag for the faint com- 
panions. The colors are directly calculated from the apparent 
magnitudes, and are not affected by parallax errors. 

The color-magnitude diagrams (with absolute magnitudes) 
for the subgroups are shown in Figure 5. The measurements 
for US are in the top panels, and those of UCL and LCC are 
in the middle and bottom panels, respectively. The curves rep- 
resent the 5 Myr (for US) and 20 Myr (for UCL and LCC) 
isochrones. The gray-shaded area indicates the inaccuracy in 
isochrone placement due to the age uncertainty in the sub- 
groups (~ 1 Myr for US members, ~ 4 Myr for UCL/LCC 
members; see Table 1). Due to our small sample and the errors 
in the photometry we cannot see a difference between the mag- 
nitude and color distributions of the three subgroups. For com- 
parison we have included the free-floating brown dwarf candi- 
dates in US reported by Martin et al. (2004). 

4.2. Distinction between companions and background 
stars 

Individual distances to Sco OB2 member stars are available (de 
Bruijne 1999), as well as the ages of the three subgroups of 
Sco OB2. A companion of a Sco OB2 member star has (by def- 
inition) practically the same distance as its primary. It is very 
likely that a primary and companion in Sco OB2 are coeval. 
The probability that a companion is captured dynamically is 
very small, since this involves either a multiple-star interac- 
tion or significant tidal dissipation. It is even less likely that the 
companion is a captured field star, i.e., that the primary and 
companion have a different age. Background stars generally 
have a different age, distance, or luminosity class than Sco OB2 
members. In principle it is possible that an ADONIS or NACO 
field contains two members of the Sco OB2 association at dif- 
ferent distances due to projection effects, but Kouwenhoven 
et al. (2005) showed that this effect can be neglected. Physical 
companion stars thus have the same age and distance as their 
primary, and therefore should fall on the isochrone for the sub- 
group to which the primary belongs, contrary to background 
stars. We use this property to separate physical companions and 
background stars. 

For each stellar component we determine in the color- 
magnitude diagram (Figure 5) the point on the isochrone that 
corresponds best to the measured position. The differences in 
color and magnitude of the star and the nearest point on the 
isochrone are then compared to the observational errors. We 
use the^- 2 test to determine how compatible the observed color 
and magnitude of a secondary are with the isochrone. For ex- 
ample, if the best-fitting value in the ( J - K s , M Ks )-diagram 
deviates by A(J-K S ) in color and by AM Ks in magnitude from 
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Fig. 4. The color-color diagram (left) and color-magnitude diagram (right) of the objects in our sample. Measurements are shown for the 
22 targets observed with NACO and for the 9 targets with multi-color observations in the ADONIS sample. Both panels show target stars 
(circles), confirmed and candidate companions (squares), and background stars (triangles). The target stars and secondaries in the US subgroup 
are indicated with filled symbols; those from UCL and LCC are indicated with open symbols. The lcr error bars are indicated for all data points. 
Lower limits are given for objects that are not detected in all three filters. Several detected objects are outside the ranges of the figures; these 
are all background stars. The status (companion or background star) of the secondaries is discussed in § 4.2. The 5 Myr isochrone for US and 
the 20 Myr isochrone for UCL and LCC are indicated with the solid and dotted curves, respectively. 



the isochrone, the x 2 value is given by 

cr, „ + cr. _ . cr, + cr 

J-K s J-K s ,iso Mk s M Ks ,iso 

where ctj-k s and <tm Ks are the observational errors in color 
and absolute magnitude, respectively. The errors on the loca- 
tion of the isochrone (due to age and metallicity uncertainty) 
are denoted with crj_ Ks is0 and cr Mjc iso . The age uncertainty is 
~ 1 Myr for members of US and ~ 4 Myr for members of 
UCL and LCC (see Table 1). The error in the placement of the 
corresponding isochrones due to age uncertainty is shown in 
Figure 5, and is small compared to the photometric errors. We 
assume a solar metallicity for all observed stellar components. 
The metallicity [M/H] of Sco OB2 has not been studied in de- 
tail. In their metallicity study of Ori OB2, Cunha & Lambert 
(1994) found a metallicity slightly (~ 0.2 dex) lower than so- 
lar for this association. Using the models of Siess et al. (2000) 
we estimate that the metallicity uncertainty A[M/H] = 0.2 re- 
sults in an additional isochronal error of 0.05 mag in JHK S and 
0.06 mag in the colors of low-mass (M 5 0.5 M Q ) companions. 

The^f 2 values for the secondaries in each of the three color- 
magnitude diagrams are listed in Table A. 2 (as long as they 
are available). For the classification into companions and back- 
ground stars we consider the largest of the three^ 2 values avail- 
able for each secondary. We choose this strategy (instead of, 
e.g., taking the average x 2 value), because background stars 
may be consistent with the isochrone for e.g. J-K s , but not for 
H - K s . The physical companions, however, should be consis- 



Secondary status Symbol Criterion 

Confirmed companion c "x 1 < 2.30 

Candidate companion ? 2.30 < x 2 — H-8 

Background star b 11.8 < x 2 

Table 3. Criteria adopted to separate the secondaries with multi-color 
observations into confirmed companions, candidate companions, and 
background stars. The^ 2 values of 2.30 and 11.8 correspond to the 
lcr and 3cr levels. This means that (statistically) 68.3% of the physical 
companions have^ 2 < 2.30 and 99.73% of the physical companions 
have^ 2 < 11.8. A secondary with^ 2 < 2.30 is very likely a compan- 
ion star. A secondary with^ 2 > 1 1.8 is almost certainly a background 
star. The secondaries with 2.30 < x 2 < H-8 may be companions or 
background stars. 



tent with the isochrone in all color-magnitude diagrams (within 
the error bars). 

Table 3 lists the criteria we adopt to classify the secondaries 
into three groups: confirmed companions, candidate compan- 
ions, and background stars. The^- 2 value of 2.30 corresponds to 
the lcr confidence level, which means that statistically 68.3% 
of the companion stars have^f 2 < 2.30. Similarly, 99.73% of 
the companions have x 2 < H-8 (3cr confidence level). The 
above confidence levels are for a dataset with two degrees of 
freedom, under the assumption that the errors are Gaussian 
(Press et al. 1992). 

For several objects we have lower limits on the magnitudes 
in one or two of the filters. Using the lower limits on J, H, or 
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Fig. 5. The (absolute) color-magnitude diagrams for the 22 targets in our NACO sample and the 9 targets with multi-color observations in the 
ADONIS sample. The results are split into the three subgroups US (top), UCL (middle), and LCC (bottom). The primary stars are indicated 
with circles; the confirmed companions with large squares, and the candidate companions with small squares. The M Ks magnitude is derived 
from the K s magnitude by correcting for distance and extinction for each target individually. The solid curves represent isochrones of 5 Myr 
(for US) and 20 Myr (for UCL and LCC). The 15 Myr and 23 Myr isochrones enclose the gray-shaded area and represent the uncertainty in 
the age of the UCL and LCC subgroups. For each data point we indicate the lcr errors. The photometry of the observed objects cannot be used 
to distinguish between the subgroups of Sco OB2, due to the errors and the small sample. The free-floating brown dwarfs in US identified by 
Martin et al. (2004) are indicated with triangles, adopting a distance of 145 pc. 



Ks we calculate upper or lower limits on x an d are able to 
classify several additional objects as background stars. 

Even though we are sensitive (although not complete) to 
massive planets around our NACO targets, we do not classify 
the very faint secondaries in this paper, for two reasons. First, 



many faint background stars are expected in the ADONIS and 
NACO field of view (see § 4.3). Due to the larger error bars for 
faint secondaries, several background stars may be consistent 
with the isochrone; the vast majority of the "candidate plan- 
ets" are likely background stars. Second, the presently available 
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K s < 12 mag 
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10 
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23 (4) 


12 < K s < 14 mag 


- (-) 


1 
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1 (1) 


K s > 14 mag 


- (-) 




(3) 
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(3) 


1 (6) 


no K s available 


- (-) 




(-) 




(1) 


- (1) 


Total 


9 (-) 


11 


(5) 


5 


(7) 


25 (12) 



Table 5. The distribution of companion stars over angular separa- 
tion p and K s magnitude for our sample of 31 targets with multi- 
color observations. Each entry lists the number of confirmed com- 
panion stars, i.e., the secondaries with x 2 < 2.30. Between brack- 
ets we list the number of candidate companion stars, which have 
2.30 < x 2 < 11-8. We have also included the candidate companion 
HIP80142-2 (p = 5.88"), for which no K s measurement is available. 
Several candidate companions are likely to be background stars, espe- 
cially faint candidates at large separation. In the region 1" < p < 4" 
and 12 mag < K s < 14 mag (which we will study in § 6) we find one 
confirmed companion and no candidate companions. 



evolutionary models for massive planets are not very reliable 
for young ages (see e.g., Chabrier et al. 2005, for a review). 
Throughout our analysis we do not consider objects with a 
mass below 0.02 M . Consequently, all secondaries with an in- 
ferred mass smaller than 0.02 M are classified as background 
stars. 

The status given to each secondary is listed in Table A. 2. 
The secondaries with^ 2 < 2.30 are very likely companions be- 
cause of their proximity to the isochrone. As statistically only 
1 out of 370 companions have x 1 > 1 1-8, we claim with high 
confidence that the secondaries with^ 2 > 1 1.8 are background 
stars. We cannot confirm the status of the secondaries with 
2.30 < x 2 < H-8. These secondaries are consistent with the 
isochrone within the 3cr error bars. However, due to the large 
number of faint background stars, it is likely that several back- 
ground stars also satisfy this criterion. As we cannot confirm 
either their companion status or their background star status, 
we will refer to the secondaries with 2.30 < x 2 < 1 1-8 as can- 
didate companions. 

The properties of the 25 confirmed companion stars found 
around the 22 members in the NACO survey and the 9 mem- 
bers with multi-color observations in the ADONIS survey are 
listed in Table 4. Table 5 lists the distribution of confirmed 
and candidate companions over angular separation p and Ks 
magnitude. The largest fraction of candidate companions (rel- 
ative to the number of confirmed companions) is seen for faint 
secondaries at large angular separation. In Section 6 we will 
study the virtual absence of companions with 1" < p < A" and 
12 < K s < 14 mag. Table 5 shows that only one confirmed 
companion is detected in this region; no candidate companions 
are found. 

4.3. The background star population 

Our method to separate companions and background stars is 
based on a comparison between the location of the secondaries 
in the color-magnitude diagrams and the isochrone. The num- 
ber of background stars identified with this method should be 
comparable to the expected number of background stars in the 



fields around the targets. In this section we make a comparison 
between these numbers, where the expected number of back- 
ground stars is based on (1) the Besancon model of the Galaxy, 
and (2) the background star study in Sco OB2 performed by 
Shatsky & Tokovinin (2002). 

We use the Besancon model of the Galaxy (Robin et al. 
2003) to characterize the statistical properties of the back- 
ground star population. We obtain star counts in the direction 
of the centers of the three subgroups, as well as for (I, b) = 
(300°, 0°), where LCC intersects the Galactic plane. We in- 
clude objects of any spectral type, luminosity class, and pop- 
ulation, up to a distance of 50 kpc, and convert the model K 
magnitude to K s magnitude. As expected, the Besancon model 
shows a strong variation in the number of background stars 
with Galactic latitude. Most background stars are found in the 
direction of the Galactic plane. For a given numerical value of 
the magnitude limit, more background stars are expected to be 
found in the Ks band than in the J and H bands. Although 
the Besancon model is in good agreement with the observed 
properties of the Galaxy, it cannot be used to make accurate 
predictions for individual lines of sight. For example, the high 
variability of the interstellar extinction with the line-of-sight, 
which is known to be important for the background star statis- 
tics (Shatsky & Tokovinin 2002), is not taken into account in 
the Besancon model. A high interstellar extinction reduces the 
observed number of background stars significantly, which is es- 
pecially important for the US subgroup, which is located near 
the p Oph star forming region. 

Let F(Ks) be the number of background stars brighter than 
K s , per unit of surface area. Above we mentioned that the num- 
ber of background stars in the Besancon model, i.e., the nor- 
malization of F(Ks), varies strongly with Galactic coordinates. 
The profile of F(K S ), however, is very similar for different lines 
of sight, and can be approximated with a function of the form 
F(K S ) = CiXl0 rKs , withy = 0.32±0.01 for5 < K s < 20 mag. 
The constant C, defines the normalization of F(K S ), which de- 
pends on the Galactic coordinates. 

The number of background stars within a certain angular 
separation p is proportional to the enclosed area A(p) in the 
field of view within that angular separation. For our NACO ob- 
servations we have a square detector of size Lnaco = 14 arc- 
sec; for the observing strategy we used for the ADONIS ob- 
servations we effectively have Lado = § • 12.76 arcsec (see 
Kouwenhoven et al. 2005). For a given p (in arcsec) the en- 
closed area (in arcsec 2 ) is then given by 



Ai(p) 



( p 2 for p < Li/2 

Tip 1 - 4p 2 arccos (L ( /2p) + L t ^4p 2 - L? 

for Li/2 <p<Li/^2 
L 2 for p > Li/ V2 



where the subscript ; refers to either the ADONIS or the NACO 
observations. Here we make the assumption that the target star 
is always in the center of the field of view. In our NACO 
survey we occasionally observe the target star off-axis in or- 
der to study a secondary at angular separation larger than 
Wo/ V2 = 9.9", but we ignore this effect here. 
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Star 


J 


H 


K s 


P 


PA 


Mj 


M H 


M Ks 


Mass 


Status 


Remarks 




mag 


mag 


mag 


arcsec 


deg 


mag 


mag 


mag 


M Q 






NACO targets 


HIP59502 -1 


12.35 


11.83 


11.64 


2.94 


26.39 


7.39 


6.86 


6.68 


0.14 


c 




HIP62026 -1 


8.08 


7.90 


7.86 


0.23 


6.34 


2.88 


2.71 


2.66 


1.19 


c 




HIP63204 -2 


8.79 


8.51 


8.40 


0.15 


236.56 


3.59 


3.31 


3.19 


1.06 


c 




HIP67260 -1 


8.88 


8.46 


8.36 


0.42 


229.46 


3.42 


2.99 


2.90 


1.10 


c 




HIP67919 -1 


9.98 


9.38 


9.10 


0.69 


296.56 


4.89 


4.30 


4.02 


0.75 


c 




HIP68532 -1 


10.52 


9.85 


9.54 


3.05 


288.50 


5.03 


4.36 


4.05 


0.73 


c 




HIP68532 -2 


11.38 


10.94 


10.63 


3.18 


291.92 


5.89 


5.45 


5.14 


0.39 


c 




HIP69113 -1 


10.98 


10.43 


10.29 


5.34 


65.15 


4.83 


4.28 


4.14 


0.77 


c 




HIP69113 -2 


11.27 


10.45 


10.30 


5.52 


67.17 


5.12 


4.29 


4.15 


0.72 


c 




HIP73937 -1 


> 8.40 


8.46 


8.37 


0.24 


190.58 


> 2.94 


3.00 


2.91 


1.11 


c 




HIP79739 -1 


12.28 


11.52 


11.23 


0.96 


118.33 


6.34 


5.58 


5.29 


0.16 


c 




HIP79771 -1 


12.00 


11.28 


10.89 


3.67 


313.38 


6.06 


5.33 


4.94 


0.19 


c 




HIP79771 -2 


12.39 


11.79 


11.42 


0.44 


128.59 


6.44 


5.85 


5.47 


0.13 


nc 




HIP80799 -1 


10.60 


10.04 


9.80 


2.94 


205.02 


5.08 


4.51 


4.27 


0.34 


c 




HIP80896 -1 


11.16 


10.63 


10.33 


2.28 


177.23 


5.60 


5.07 


4.77 


0.24 


c 




HIP81972 -3 


12.54 


11.86 


11.77 


5.04 


213.45 


6.16 


5.48 


5.39 


0.35 


c 


J 


HIP81972 -4 


15.10 


14.43 


13.98 


2.79 


106.94 


8.72 


8.05 


7.60 


0.06 


nc 


JHK 


HIP81972 -5 


16.11 


15.63 


15.26 


7.92 


229.27 


9.73 


9.25 


8.88 


-0.03 


nc 


JHK 


ADONIS multi-color subset 


HIP76071 -1 


> 11.25 


11.28 


10.87 


0.69 


40.85 


> 5.09 


5.12 


4.71 


0.23 


c 




HIP77911 -1 


12.68 


12.20 


11.84 


7.96 


279.25 


6.82 


6.34 


5.98 


0.09 


c 




HIP78809 -1 


11.08 


10.45 


10.26 


1.18 


25.67 


5.32 


4.69 


4.50 


0.30 


c 




HIP78956 -1 


9.76 


9.12 


9.04 


1.02 


48.67 


3.39 


2.75 


2.67 


1.16 


c 




HIP79124 -1 


11.38 


10.55 


10.38 


1.02 


96.18 


5.33 


4.50 


4.33 


0.33 


c 




HIP79156 -1 


11.62 


10.89 


10.77 


0.89 


58.88 


5.50 


4.77 


4.65 


0.27 


c 




HIP80238 -1 


7.96 


7.66 


7.49 


1.03 


318.46 


2.34 


2.04 


1.87 


1.67 


c 





Table 4. Properties of the 25 confirmed companion stars found around the 22 members in our NACO survey and the 9 members with multi- 
color observations in the ADONIS survey. The columns show the secondary designation, the J, H, and K s magnitudes, the angular separation, 
and the position angle (measured from North to East). Magnitude lower limits are given if a secondary is not detected in a filter. We list the 
absolute magnitude and corresponding mass in columns 7 - 10. The 11th column lists the status of the object (c = confirmed companion star, 
nc = new confirmed companion star). The last column shows additional remarks. A "J", "H", or "K" means that the secondary flux in this filter 
was obtained from the image obtained without the NDF, using the PSF from the corresponding image that was obtained with NDF (see § 2.4). 
Properties of the observed primaries, candidate companions, and background stars are not shown here; these are listed in Table A.l. 



We now have expressions for the quantity N(K s ,p), the ex- 
pected number of background stars with magnitude brighter 
than K s and angular separation smaller than p, as a function 
of K s and p: 

Ni(K s ,p) = F(K S )-Ai(p) = Q ■ W Ks ■ Ai(p) , (3) 

where y = 0.32 + 0.01, A(p) is defined in Equation 2, and C, is 
a normalization constant. 

For the normalization of Equation 3 we compare our obser- 
vations with the background star study of Shatsky & Tokovinin 
(2002). Apart from their target observations, the authors addi- 
tionally obtained sky images centered at 21 arcsec from each 
target in order to characterize the background star popula- 
tion. From the background star study of Shatsky & Tokovinin 
(2002) we derive the expected number of background stars in 
the ADONIS and NACO field of view. We assume a detection 
limit of K s = 15 mag for ADONIS and K s = 17 mag for 
NACO (which roughly corresponds to the completeness limit 
of the background star study of Shatsky & Tokovinin 2002). 
Using the data from the background star study of Shatsky & 
Tokovinin (2002) we calculate the expected number of back- 
ground stars for each ADONIS and NACO field, for the three 



ADONIS NACO 



Region 


stars per field 


^fields 


stars per field 


Affields 


US 

UCL 

LCC 


06 + " 06 
u.uu „ 02 

f> 97 +0.06 
KJ - Z - 1 -0.02 

0.43 ± 0.08 


51 
64 
40 


10 +0 11 
47 +<m 
0.73 ±0.14 


9 

3 
5 


GP 


1.51 ± 1.03 


22 


2.59 ± 1.76 


5 


Total 




177 




22 



Table 6. The number of background stars expected in our ADONIS 
and NACO fields of view, based on the background star study of 
Shatsky & Tokovinin (2002). Column 1 lists the four regions for which 
we study the background statistics. The targets with Galactic latitude 
\b\ < 5° are included in the region GP. The other targets are grouped 
according to their membership of US, UCL, and LCC. For the 177 tar- 
gets only observed with ADONIS and the 22 targets observed with 
NACO we list the expected number of background stars per field of 
view, and the number of targets observed in the four regions. In to- 
tal we expect to find 70.75 ± 4.90 background stars in the ADONIS 
sample, and 18.96 ± 3.96 in the NACO sample. 



subgroups and for targets close to the Galactic plane Qb\ < 5°). 
Table 6 lists the expected number of background stars with cor- 
responding Poisson errors. Additionally, we list the number of 



Kouwenhoven et al.: A brown dwarf desert for intermediate mass stars in Sco OB2? 



13 




Fig. 6. The expected number of background stars brighter than K s at 
angular separation smaller than p for the combined sample of 199 tar- 
gets (solid contours). Overplotted are the background stars detected 
around the 22 NACO targets (filled triangles), and those found around 
the 177 targets observed with ADONIS only (open triangles). The ob- 
served cumulative number of background stars is indicated with the 
dotted contours, for values of 10 and 50 background stars, respec- 
tively. In Kouwenhoven et al. (2005) we classify secondaries brighter 
than Ks = 12 mag (upper horizontal line) as candidate companions, 
and those fainter than K s = 12 mag as probable background stars. The 
faintest objects we detect in our ADONIS and NACO surveys have 
Ks magnitude of approximately 15 and 17 magnitudes, respectively 
(lower horizontal lines). 



targets observed in each of the four regions. For the 177 tar- 
gets observed with ADONIS only we expect « 71 background 
stars, and for the 22 targets observed with NACO we expect 
« 19 background stars. In our combined ADONIS and NACO 
dataset, the expected number of background stars is 90+6. This 
gives normalization factors Cadonis = 1.74xl0~ 8 arcsec~ 2 and 
Cnaco = 1.60 x 10~ 8 arcsec~ 2 for Equation 3. 

Figure 6 shows the expected number of background stars 
brighter than Ks and closer than p as a function of K s and 
p, for the combined ADONIS and NACO dataset. Since most 
of our targets are observed with ADONIS only, the shape 
of the expected background star density distribution is domi- 
nated by that of ADONIS. As we have a square field of size 
Lado = 19.1 arcsec, the cumulative number of background 
stars rises rapidly between p = 0" and p = Lado/ V2 = 13.5" 
and becomes flat for larger separations. The background stars 
from the NACO survey are represented with the filled triangles, 
and those detected only in the ADONIS survey are represented 
with open triangles. Figure 6 shows that the (p,K s ) distribu- 
tion of the observed 97 background stars is in good agreement 
with that of the expected 90 background stars. In the extreme 
case that all 12 candidate companions are actually background 



Fig. 7. The x 2 distance to the isochrone in the (/ - K s , M Ks ) diagram 
versus N(K s ,p), where N(K s ,p) is the expected number of back- 
ground stars brighter than Ks and angular separation smaller than 
p. The vertical dashed lines correspond to x 2 = 2.30 (lcr; left) and 
X 2 = 11.8 (3<x; right). The symbols represent the confirmed com- 
panions (filled squares), the candidate companions (open squares), 
and background stars (triangles). Note that the classification of a sec- 
ondary is based on the x 2 values for the different color-magnitude 
diagrams; not only for the x 2 of the (J - K s , M Ks )-diagram, which 
is shown above. The expected number of background stars in a field 
(vertical axis) is used as a consistency check. The horizontal dotted 
line represents the 1% filter used by Poveda et al. (1982) to sepa- 
rate companion stars (below the line) and background stars (above 
the line). The 1% filter is a reasonable method when multi-color ob- 
servations are not available, but is not used in our study. The triangle 
in the lower-right quadrant represents HIP63204-1 (see § 4.4). 



stars, the observed number of background stars is 110. In this 
extreme case there are 22% more background stars than ex- 
pected, which suggests that a significant part of the candidate 
companions may indeed be physical companions. 

Background stars are by definition not associated with the 
target star, and therefore generally have (1) a random position 
in the field of view of the observation, and (2) a position in the 
color-magnitude diagram that is likely to be inconsistent with 
the isochrone. Figure 7 shows the relation between N(Ks,p) 
and the x 1 value derived from the location of the secondary in 
the ( J - K s , M Ks ) diagram with the isochrone . In other words, 
Figure 7 shows the probability of detecting a background star at 
separation p (or smaller) and magnitude Ks (or brighter) versus 
how far away the secondary is from the isochrone. The verti- 
cal dashed lines in Figure 7 are at^ 2 = 2.30 and^ 2 = 11.8, 
the values used to determine the status of the companions (see 
Table 3). This correlation provides additional support to the 
method we use to separate companions and background stars. 
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Poveda et al. (1982) performed a statistical study of bi- 
nary stars in the Index Catalogue of Visual Double Stars. They 
showed that it is statistically plausible to assume that compo- 
nents with N(p,rri2) > 0.01 are background stars, where ni2 is 
the magnitude of the secondary. This technique is referred to as 
the "1% filter". The horizontal line in Figure 7 represents the 
1% filter used by Poveda et al. (1982). Secondaries below this 
line would be classified as companions using the 1 % filter, and 
those above would be classified as background stars. Figure 7 
shows that the 1 % filter is a reasonable technique, but not as 
accurate as the multi-color technique used in this paper. 

4.4. Notes on some individual secondaries 

We detect two hierarchical triple systems: HIP68532 and 
HIP69113. The two companions of HIP68532 (previously re- 
ported in Kouwenhoven et al. 2005) have a mass of 0.73 M 
and 0.39 M Q , respectively. HIP68532 has a companions-to- 
primary mass ratio of (0.73 + 0.39)/ 1.95 = 0.57. The compan- 
ions are separated 0.23" (~ 28 AU) from each other and 3.11" 
(~ 385 AU) from the primary, giving an estimate of 0.073 for 
the ratio between the semi-major axes of the inner and outer 
orbits. The two companions of HIP691 13 (previously reported 
in Huelamo et al. 2001) have a mass of 0.77 M and 0.72 M , 
respectively, corresponding to a companions-to-primary mass 
ratio of 0.39. The companions are separated 0.26" (~ 44 AU) 
from each other and 5.43" (~ 917 AU) from the primary, giv- 
ing an estimate of 0.048 for the ratio between the semi-major 
axes of the inner and outer orbits. 

For HIP62026-1 we find a significant difference in po- 
sition angle between the ADONIS and NACO observations. 
With the ADONIS observations, obtained on 8 June 2001, we 
find (p,<p) = (0.22", 12.5°). With NACO we measure (p,<p) = 
(0.23", 6.34°), 1033 days later. As the angular separation be- 
tween HIP62026-1 and its primary is small, the observed po- 
sition angle difference may well be the result of orbital mo- 
tion. Assuming a circular, face-on orbit, we estimate an or- 
bital period of 165 year for the system HIP62026. The total 
mass of the system (taken from Table A.l) is 3.64 ± 0.25 M Q , 
which gives via Kepler's law a semi-major axis of 46 AU. This 
value is of the same order of magnitude as the (projected) semi- 
major axis of ~ 24 AU derived from the angular separation of 
0.22" between the components (adopting a distance of 109 pc 
toHIP62026). 

HIP63204-1 is a bright and red object separated only 
1.87 arcsec from the LCC member HIP63204. The isolated lo- 
cation of HIP63204-1 in the bottom-right quadrant of Figure 7 
shows that the probability of finding a background star of this 
magnitude (or brighter) at this angular separation (or closer) 
is small. According to its location in the color-magnitude di- 
agrams, HIP63204-1 is a background star and hence we clas- 
sify it as such. HIP63204 and its companion HIP63204-2 at 
p = 0.15" have masses of 2.05 M Q and 1.06 M Q , respectively. 
If HIP63204-1 (at p = 1.87") would be a companion, its mass 
would be approximately 1 M , in which case HIP63204 would 
be an unstable triple system. The colors of HIP63204-1 are con- 
sistent with a 0.075 M Q brown dwarf with an age of 10 Gyr at 



a distance of 60 pc, and are also consistent with those of an 
M5 III giant at a distance of ~ 5.6 kpc (using the models of 
Cox 2000; Chabrier et al. 2000). 

HIP8 1972-3, HIP8 1972-4, and HIP8 1972-5 fall on the 
20 Myr isochrone in all three color-magnitude diagrams. These 
objects are likely low-mass companions of HIP81972 (see 
§ 5). HIP8 1972-5 is the topmost companion (black square) in 
Figure 7. HIP8 1972-5 is the faintest companion in our sample, 
and the expected number of background stars with a similar or 
brighter magnitude and a similar or smaller separation is large 
(~ 16 for the ADONIS sample). The secondary HIP8 1972-2 
(at p = 7.02") was reported before as a possible companion of 
HIP81972 in Worley & Douglass (1997), but the secondary is 
classified as a background star by Shatsky & Tokovinin (2002). 
With our NACO multi-color observations we cannot determine 
the nature of this secondary with high confidence. The LCC 
member HIP81972 is located close to the Galactic equator 
(b = +3° 11'), so care should be taken; background stars with 
a magnitude similar to that of the secondaries of HIP8 1972 are 
expected in the field around this star. 

Kouwenhoven et al. (2005) identified seven "close back- 
ground stars" (Ks > 14 mag; p < 4"), for which the back- 
ground star status was derived using the K s magnitude only. 
These are objects next to the targets HIP61265, HIP67260, 
HIP73937, HIP78958, HIP79098, HIP79410, and HIP81949. 
The ADONIS observations are incomplete in this region (see 
Figure 3). More low-mass companions with Ks > 14 mag and 
1" < p < 4" may be present for the 177 A and late-B stars only 
observed with ADONIS. 

With our NACO multi-color observations we confirm 
that five of the "close background stars" (HIP61265-2, 
HIP73937-2, HIP79098-1, 79410-1, and HIP81949-2) are 
background stars. For the other two secondaries, HIP67260-3 
and HIP78969-1, we cannot determine whether they are back- 
ground stars or brown dwarf companions. As many background 
stars with similar magnitudes are expected in the field, these are 
likely background stars. However, follow-up spectroscopic ob- 
servations are necessary to determine the true nature of these 
close secondaries. 

4.5. Accuracy of the K s = 12 separation criterion 

One of the goals of our study is to evaluate the accuracy of 
the Ks = 12 mag criterion that we used to separate compan- 
ions and background stars in the ADONIS survey. This is pos- 
sible, now that we have performed a multi-color analysis of 
72 secondaries around 31 members of Sco OB2. Table 7 shows 
the distribution of secondary status for the secondaries with 
K s < 12 mag and those with K s > 12 mag (three secondaries 
without K s measurements are not included). According to the 
Ks = 12 criterion, all secondaries brighter than Ks = 12 mag 
are companions, while all fainter secondaries are background 
stars. 

If we consider only the confirmed companions and back- 
ground stars, we see that the Ks = 12 mag criterion correctly 
classifies the secondaries in / = (23 + 27)/(25 + 33) = 86% 
of the cases. If all candidate companions are indeed compan- 
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Status 


K s 


< 12 mag 


K s 


> 12 mag 


Total 


c 


23 


(70%) 


2 


(6%) 


25 


? 


4 


(12%) 


7 


(19%) 


11 


b 


6 


(18%) 


27 


(75%) 


33 


Total 


33 


(100%) 


36 


(100%) 


69 



Table 7. Accuracy of the Ks = 12 mag criterion to separate compan- 
ions and background stars. This table contains 69 out of the 72 secon- 
daries with multi-color observations in the ADONIS or NACO dataset. 
Three secondaries (1 candidate companion; 2 background stars) for 
which no K s magnitudes are available, are not included. The first col- 
umn shows the status of the secondary (c = confirmed companion, ? 
= candidate companion, b = background star). Columns 2 to 5 list 
the distribution over status for secondaries with K s < 12 mag and 
K s > 12 mag. Depending on the true nature of the candidate com- 
panion stars, the K s = 12 mag criterion correctly classifies the secon- 
daries in ~ 80% of the cases. 



ions, we have / = 78%, while if all candidate companions are 
background stars, we have / = 82%. This indicates that ~ 80% 
of the candidate companions identified by Kouwenhoven et al. 
(2005) are indeed companion stars. 

The K s = 12 mag criterion is accurate for the measured 
set of secondaries as a whole. It is obvious that for the lowest- 
mass companions the criterion is not applicable, as virtually all 
brown dwarf and planetary companions have Ks > 12 mag at 
the distance of Sco OB2. Out of the 25 confirmed companions 
we find with our multi-color analysis, 23 indeed have K s < 
12 mag, but two have K s > 12 mag. These are the brown dwarf 
companions of HIP81972 (see§ 4.2 and 5 for a discussion). 

Now that we have confirmed the validity of the K s = 
12 mag criterion, many of the candidate companions found in 
the ADONIS survey very likely are physical companion stars. 
Table A. 3 gives an overview of all candidate and confirmed 
companion stars identified in the ADONIS and NACO surveys. 

5. Masses and mass ratios 

For each primary and companion star we derive the mass us- 
ing its color and magnitude. We find the best-fitting mass by 
minimizing the x 1 difference between the isochrone and the 
measurements, while taking into account the errors in the mea- 
surements: 



2 (KJ-Ks) 
X = 

0~J-K s 



A(H-K S ) 

C H-K< 



(4) 



The masses of all primaries and confirmed companions are 
listed in Table A. 1 . We additionally list the masses of the can- 
didate companions, assuming that they are indeed companions, 
but we do not include these masses in our analysis. We find 
primary masses between 1.1 M Q and 4.9 M . The confirmed 
companion star masses range between 0.03 M and 1.19 M G , 
with mass ratios 0.006 < q < 0.55. The average error in the 
mass as a result of the error in the color and magnitude is 8.5% 
for the primaries and 12% for the companion stars. The aver- 
age error in the mass ratio is 12.5%. Although accurate B and 
B-V measurements are available for the primaries, we do not 
use these. The B and V measurements often include the flux 



of unresolved close companions, and therefore lead to overes- 
timating the primary masses. 

Kouwenhoven et al. (2005) derived masses from K s mag- 
nitudes only. For the primary stars these are close to those ob- 
tained from multi-color observations in our current analysis. 
The rms difference between the masses derived using the two 
methods is 6.6%. No systematic difference is present for the 
primaries. Kouwenhoven et al. (2005) overestimated the com- 
panion star masses with ~ 2.2% (~ 0.01 M ) on average. 

The companions with the lowest mass are those of 
HIP81972, which have masses of 0.35 M (370 Mj), 0.06 M 
(63 Mj), and 0.03 M (32 Mj). The latter two are likely brown 
dwarfs. With an angular separation of p = 2.79", the 63 Mj 
component is the only observed brown dwarf in the 1 " - 4" an- 
gular separation interval (see § 6). The other two companions 
of HIP81972 have a larger angular separation. 



6. The lower end of the companion mass 
distribution 

Kouwenhoven et al. (2005) discussed the potential lack of sub- 
stellar companions to A and late-B members of Sco OB2. With 
our NACO follow-up observations we confirm the very low 
number of brown dwarf companions with respect to the num- 
ber of stellar companions found around these stars (§ 6.1). In 
Section 6.2 we will discuss whether a brown dwarf desert exists 
among A and late-B members of Sco OB2. In Section 6.3 we 
will briefly discuss the potential origin of such a brown dwarf 
desert. We will show that the small brown dwarf companion 
fraction among A- and B stars in Sco OB2 can be explained 
by an extrapolation of the stellar companion mass distribution, 
i.e., there is no need to eject brown dwarf companions from 
binary systems at an early stage of the formation process (the 
embryo ejection scenario; Reipurth & Clarke 2001). 

6.1. The brown dwarf "gap" for 1" < p < 4" 

Kouwenhoven et al. (2005) observed a gap in the (p, K$) distri- 
bution of the stellar companions in the Sco OB2 binary popula- 
tion: no secondaries with a magnitude 12 mag < Ks < 14 mag 
and an angular separations p < 4" were detected. These sec- 
ondaries should have been detected, had they existed, since 
the ADONIS survey is almost complete in this region (see 
Figure 3). Figure 8 shows the distribution of Ks and AKs as 
a function of p for the ADONIS and NACO observations com- 
bined. The "gap" in the (p, K s ) distribution described above is 
clearly visible. With our NACO survey we detect one compan- 
ion at the bottom of this region: the brown dwarf companion 
HIP8 1 972-4 (K s = 13.98 ±0.12 mag; see also § 4.2 and 5). No 
other secondaries are present in this region. 

The stellar companion fraction is the fraction of stars with 
stellar companions. Among A and late-B stars in Sco OB2 in 
the semi-major axis range l"-4" (130-520 AU) we find a stel- 
lar companion fraction of 14 ± 3%. Similarly, the brown dwarf 
companion fraction is the fraction of stars with brown dwarf 
companions. The brown dwarf companion fraction for the stars 
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Fig. 8. Companion star magnitude K s (top) and magnitude difference 
AK S (bottom) versus angular separation for the combined ADONIS 
and NACO datasets. The symbols represent confirmed companions 
(filled squares), candidate companions (open squares), and back- 
ground stars (triangles); see Section 4.2 for further information on 
this classification. The horizontal line corresponds to K s = 12 mag, 
the criterion used by Kouwenhoven et al. (2005) to separate compan- 
ion stars and background stars. Typical detection and completeness 
limits corresponding to the observations are shown in Figure 3. For a 
given K s magnitude, the number of background stars closer than an- 
gular separation p is given by Equation 3. This figure clearly shows the 
dearth of brown dwarf companions with 1" < p < 4" around A and 
late-B stars in Sco OB2 (region indicated with the dashed rectangle). 



in this separation range is 0.5 ± 0.5% 4 (for brown dwarfs with 



4 Note that we only find one brown dwarf companion with 
12 mag < K s < 14 mag. We find no background stars in this region, 
so accidental misclassification of companions as background stars is 
not an issue here. 



K s < 14 mag). The substellar-to-stellar companion ratio R is 
defined as 



R 



number of brown dwarf companions 
number of stellar companions 



(5) 



In our study we cannot calculate R because we do not know 
how many faint (K s > 14 mag) brown dwarf companions 
are missing. We therefore calculate the restricted substellar-to- 
stellar companion ratio using only the brown dwarf com- 
panions brighter than K$ = 14 mag 5 . In the angular separation 
range 1"— 4" we observe R+jm = 0.036±0.036 for intermediate 
mass stars in Sco OB2. We cannot make similar statements for 
companions with properties other than those described above. 
For angular separations smaller than 1" our survey is signifi- 
cantly incomplete for 12 mag < Ks < 14 mag. For p > 4" 
many objects with 12 mag < K s < 14 mag are likely back- 
ground stars, of which the status still needs to be confirmed. 
Finally, the ADONIS survey is incomplete for K s > 14 mag. 

6.2. A real brown dwarf desert? 

The brown dwarf desert is defined as a deficit (not necessarily 
a total absence) of brown dwarf companions, either relative to 
the frequency of companion stars or relative to the frequency 
of planetary companions (McCarthy & Zuckerman 2004). In 
this paper the brown dwarf desert for A and late-B members of 
Sco OB2 is characterized by a small number of observed com- 
panions Af*,BD with 12 mag < K s < 14 mag and I" < p < 4" 
and a small substellar-to-stellar companion ratio (where the 
star in the subscript refers to the observed quantities). In gen- 
eral, the quantities A^bd and R depend on (1) the mass distri- 
bution, (2) the pairing properties of the binary systems, and 
(3) the spectral type of the stars in the sample. Unlike R, the 
value of A^bd also depends on (4) the multiplicity fraction F M 
(Equation 9), and (5) the semi-major axis (or period) distri- 
bution. We use simulated observations and compare the ob- 
served values with those predicted for various models, in order 
to roughly estimate the mass distribution and pairing proper- 
ties. For comparison between observations and simulations we 
only consider the companions brighter than K s = 14 mag in 
the angular separation range 1" < p < 4". 

We simulate models using the STARLAB package (e.g., 
Portegies Zwart et al. 2001). The primary mass is drawn in the 
mass range 0.02 M < M < 20 M , either the Salpeter mass 
distribution 

fmiM) = ^ c* M- 235 , (6) 
or from the the extended Preibisch mass distribution 



dM 

/m(M) = — «< 



M a 


for 


0.02 <M/M Q 


<0.08 


M~ 


0.9 


for 


0.08 < M/M 


<0.6 


M~ 


2.8 


for 


0.6 < M/M 


<2 


M~ 


2.6 


for 


2 < M/M 


<20 



(7) 



5 In this section we denote the observed quantities with a star as a 
subscript. For example, R denotes the substellar-to-stellar companion 
ratio (including all brown dwarfs), while R t indicates the observed 
substellar-to-stellar companion ratio, including only the brown dwarfs 
brighter than K s = 14 mag. 
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# 


Model 




N *,BD,IM 


R 


»,IM 


^V*,BD,LM 


^*,LM 


U 


ADONIS/NACO observations 


1 ± 


1 


0.036 


± 0.036 


unknown 


unknown 


1 


extended Preibisch MF, a = 


-0.9, random pairing 


5.50 ± 


0.48 


0.34 


+ 0.03 


7.19 ±0.17 


0.50 ± 0.01 


2 


extended Preibisch MF, a = 


-0.3, random pairing 


4.50 ± 


0.41 


0.24 


±0.03 


5.08 ±0.13 


0.30 ± 0.01 


3 


extended Preibisch MF, a = 


+2.5, random pairing 


1.07 ± 


0.18 


0.05 


±0.01 


1.42 ±0.07 


0.07 ± 0.01 


4 


Salpeter MF, random pairing 


r 


15.31 ± 2.79 


6.00 ± 2.90 


17.18 ±0.88 


3.95 ± 0.45 


5 


extended Preibisch MF, a = 


-0.9, f q (q) oc q-° 3i 


0.72 ± 


0.24 


0.04 


±0.01 


3.42 ±0.14 


0.18 ±0.01 


6 


extended Preibisch MF, a = 


-0.3, f q (q) oc 9 "0« 


0.71 ± 


0.22 


0.04 


±0.01 


3.35 ±0.14 


0.18 ±0.01 


7 


extended Preibisch MF, a = 


+2.5, f q (q)ocq-^ 


1.19 + 


0.27 


0.06 


±0.01 


3.30 ±0.13 


0.18 ±0.01 


8 


Salpeter MF, f q (q) oc q-°- 3i 


1.00 ± 


0.57 


0.05 


±0.02 


3.70 ±0.57 


0.20 + 0.03 



Table 8. The observed and expected number of brown dwarfs with 1" < p < 4" and 12 < K s < 14 mag for the sample of 199 target stars. 
The left columns shows the various models for which we simulated observations. Each model has a semi-major axis distribution f a (a) oc a' 1 
with 15 R < a < 5 x 10 6 R Q , and a multiplicity fraction of F M = 100%. Columns 3 and 4 show for a survey of intermediate mass stars 
(late-B and A stars; 1.4 M Q < M < 7.7 M Q ) the expected number of brown dwarfs A^.bdjm and the substellar-to-stellar companion ratio i?*,iM, 
both with l<x errors. By comparing the predicted values of A^bdjm ar, d ^*,im with the observations, we can exclude models 1, 2, and 4. In 
Kouwenhoven et al. (2005) we exclude random pairing from the Preibisch mass distribution (models 1-3) since these models are inconsistent 
with the observed mass ratio distribution of stellar companions. We additionally list the values A^,bd,lm and ^*,lm that are expected for a survey 
amongst 199 low-mass stars (0.3 M Q < M < 1.4 M Q ) in columns 5 and 6. For models with F M < 100% the expected number of brown dwarfs 
reduces to Fm x A^.bd, while R remains unchanged. Models with a smaller semi-major axis range and models with the log-normal period 
distribution found by Duquennoy & Mayor (1991) have a larger expected value of A^bd.im, A'* > bd,lm. Under the assumption that companion 
mass and semi-major axis are uncorrelated, the values of R*jm and ^*_ L m are equal to those listed above, if the log-normal period distribution 
is chosen. 



The extended Preibisch mass distribution (see Kouwenhoven 
et al. 2005) is based on the study by Preibisch et al. (2002), 
who derived /m(M) with M > 0.1 M Q for the US subgroup of 
Sco OB2. Since our current knowledge about the brown dwarf 
population in OB associations (particularly Sco OB2) is incom- 
plete (e.g., Table 2 in Preibisch et al. 2003) we simulate asso- 
ciations with three different slopes for the mass distribution in 
the brown dwarf regime. We extend the Preibisch mass distri- 
bution down to 0.02 M with a = -0.9, a = -0.3, or a = +2.5. 
The mass distribution with a = -0.9 has the same slope in 
the brown dwarf regime as for the low-mass stars. The simu- 
lations with a = -0.3 and a = +2.5 bracket the values for a 
that are observed in various clusters and the field star popula- 
tion (see Preibisch et al. 2003, for a summary). The companion 
mass is obtained via randomly pairing the binary components 
from the mass distribution or via a mass ratio distribution of 
the form f q {q) oc q~ 033 with < q < 1 and the requirement 
that any companion has a mass larger than 0.02 M . The lat- 
ter mass ratio distribution was derived from the observed mass 
ratio distribution in our ADONIS survey (Kouwenhoven et al. 
2005). For the models with random pairing, the primary star 
and companion mass are drawn independently from the mass 
distribution, and switched, if necessary, so that the primary is 
the most massive star. 

Each simulated association consists of 100000 binaries, 
has a distance of 130 pc and a homogeneous density distri- 
bution with a radius of 20 pc, properties similar to those of 
the subgroups in Sco OB2. We assume a thermal eccentric- 
ity distribution and a semi-major axis distribution of the form 
fa(a) oc a~ l , which is equivalent to /j ogfl (loga) = constant 
(Opiks law). The lower limit of a is set to 15 R . The upper 
limit is set to 5 x 10 6 R ~ 0.1 pc, the separation of the widest 
observed binaries in the Galactic disk (e.g., Close et al. 1990; 
Chaname & Gould 2004). 



Table 8 lists for eight models the predicted value of A^* jBD , 
the expected number of brown dwarfs with 1 " < p < 4" and 
12 < K s < 14 mag, normalized to a sample of 199 stars. The 
table also lists R+, the ratio between the number of brown dwarf 
companions with Ks < 14 mag and the number of stellar com- 
panions in the separation range 1" - 4". Table 9 lists A^d (the 
intrinsic number of brown dwarfs with 1" < p < 4"), and 
corresponding ratio R. In this table all companions are taken 
into account, including those with K s > 14 mag. The values in 
Table 8 can be compared directly with the observations, while 
those in Table 9 represent the intrinsic properties of each as- 
sociation model. A brown dwarf with K s = 14 mag in the US 
subgroup has a mass of less than 0.02 M (21 Mj), and a brown 
dwarf with a similar brightness in the UCL and LCC subgroup 
has a mass of ~ 0.038 M (40 Mj). 

The definition of the brown dwarf desert given in the be- 
ginning of this section is generally used for binarity studies of 
late-type stars. In our study the primaries are intermediate mass 
stars, allowing companion stars over a larger mass range than 
for low mass primaries. This naturally leads to lower values 
of A^bd and R. We therefore also list the results for a simu- 
lated survey of 199 low mass stars. Tables 8 and 9 show that 
indeed the expected values A^bd and R for low-mass stars are 
higher than those for intermediate-mass stars by ~ 30% for the 
random pairing models, and by ~ 250% for the models with 

A multiplicity fraction of Fm = 100% is assumed in each 
model. For A and late B members of Sco OB2 the observed 
multiplicity fraction F M is « 50% (Kouwenhoven et al. 2005). 
This is a lower limit of the true multiplicity fraction due to the 
presence of unresolved companions, and hence we have 50% ^ 
F M < 100%. For a multiplicity fraction smaller than 100%, the 
expected number of brown dwarfs is given by .FmxA^bd, while 
the values of R remain unchanged. In each model we adopted 
Opiks law, with 15 R < a < 5 x 10 6 R . In reality, the upper 
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# 


A'bd.im 




A'bd.lm 







unknown 


unknown 


unknown 


unknown 


1 


6.68 ± 0.53 


0.44 ± 0.04 


8.50 ±0.18 


0.65 ± 0.02 


2 


5.46 ± 0.45 


0.31 ±0.03 


6.42 ±0.15 


0.42 ± 0.01 


3 


2.01 ± 0.25 


0.10 ±0.01 


2.65 ± 0.09 


0.14 ±0.01 


4 


15.82 ±2.84 


7.75 ±4.12 


18.26 ±0.91 


5.60 ± 0.72 


5 


1.20 ± 0.31 


0.06 ± 0.02 


4.26 ±0.16 


0.24 ± 0.01 


6 


1.13 ± 0.28 


0.06 ± 0.02 


4.30 ±0.16 


0.25 ± 0.01 


7 


1.42 + 0.29 


0.07 ± 0.01 


4.18 ±0.14 


0.24 ±0.01 


8 


1.19 ±0.63 


0.06 ± 0.03 


5.02 ± 0.66 


0.30 ± 0.04 



Table 9. The observed and expected number of brown dwarfs with 
1" < p < 4" for the sample of 199 target stars. In this table we show 
the results for the full brown dwarf mass range (0.02 M Q < M < 
0.08 M Q ), unlike in Table 8, where we show the results for the brown 
dwarfs restricted to 12 < K s < 14 mag. 



limit for a may be smaller, as Sco OB2 is an expanding as- 
sociation (Blaauw 1964; Brown et al. 1999). If this is true, the 
values for 7V*3d are underpredicted, as less companions are ex- 
pected to have very large separations. Furthermore, instead of 
Opiks law, it may also be possible that the log-normal period 
distribution found by Duquennoy & Mayor (1991) holds. For a 
model with Opiks law at a distance of 130 pc, 11% of the com- 
panions have separations between 1 " - 4", while for a model 
with the log-normal period distribution, 13% of the compan- 
ions have separations between 1" - 4". If the log-normal pe- 
riod distribution holds, the values for /V* >BD in Tables 8 and 9 
are underpredicted. The value of R (and /?*) does not change 
for the possible adaptations described here, under the assump- 
tion that the stellar and substellar companions have the same 
semi-major axis (or period) distribution. 

Table 8 shows that the models with f q {q) °c q~ 033 are in 
good agreement with our observations for any value of a. The 
reason for this is that A^bd.im and R+jm are independent of a 
for these models, as only the primary is chosen from the mass 
distribution. For the models with random pairing, the two com- 
ponents of each binary system are independently chosen from 
the mass distribution. Only those models with a turnover in the 
mass distribution in the brown dwarf regime are consistent with 
the observations (for a multiplicity fraction of 0.5 $ Fm < 1). 
However, in Kouwenhoven et al. (2005) we excluded random 
pairing by studying the observed mass ratio distribution for 
stellar companions. The remaining models that are consistent 
with our observations have an extended Preibisch mass distri- 
bution and a mass ratio distribution of the form f q {q) oc g~ 033 . 
Although this distribution is peaked to low values of q, the 
number of brown dwarf companions is significantly smaller 
than the number of stellar companions. For example, for a sam- 
ple of binaries with a primary mass of 3 M , the substellar-to- 
stellar companion mass ratio R (see Equation 5) resulting from 
fq(q) K q~ 033 is given by 

X 0.08/3 
f q (q)dq r 0.67 1 °.°8/3 
^ ,2/3 = [g J 0.02/3 =0Q59 

p r «°- 67 1 

fq(q)dq Y q J 0.08/3 

Jo.08/3 



Figure 9 further illustrates that a small value for R is expected 
among binaries with an intermediate-mass or solar-type pri- 
mary, even if f q (q) is strongly peaked to low values of q. 
Among primaries with a mass of 1 M and 3 M , this fraction is 
~ 14% and ~ 6%, respectively. For this mass ratio distribution, 
the number of brown dwarf companions is significantly smaller 
than the number of stellar companions, even if observational bi- 
ases are not taken into account. If binary formation truly results 
in a mass ratio distribution similar to f q (q) oc q~ 033 , the brown 
dwarf desert (in terms of the "deficit" of brown dwarf compan- 
ions relative to stellar companions) is a natural outcome of the 
star forming process for intermediate mass stars. 

The observed number of brown dwarfs (with 12 mag < 
Ks < 14 mag) is A^bdjm = 1 ± 1. After correction for unseen 
low-mass brown dwarfs (with Ks > 14 mag) this translates to 
A^bd.im = 1.6 ± 1.6 brown dwarfs (cf. Tables 8 and 9). If we 
assume a semi-major axis distribution of the form f a (a) °c cT l 
with 15 R© < a < 5 x 10 6 R Q and a distance of 130 pc, we 
expect ~ 1 1 % of the companions to be in the angular separa- 
tion range 1" - 4". Assuming that companion mass and semi- 
major axis are uncorrelated, this also means that 11% of the 
brown dwarfs are in this range. Extrapolation gives an estimate 
of (1.6 ± 1.6)/0.11 = 14.5 + 14.5 brown dwarf companions 
around the 199 target stars, or a brown dwarf companion frac- 
tion of 7.3 ± 7.3% for intermediate mass stars in Sco OB2. On 
the other hand, if we assume the log-normal period distribution 
found by Duquennoy & Mayor (1991), we find a correspond- 
ing brown dwarf companion fraction of 6.2 ±6.2% (for a model 
binary fraction of 100%). Note that if a mass ratio distribution 
f q (q) is adopted, these values are independent of the slope a of 
the mass distribution in the brown dwarf regime. 

McCarthy & Zuckerman (2004) find a companion fre- 
quency of 0.7±0.7% for brown dwarf companions (M > 30 Mj) 
to F, G, K, and M stars in the separation range 120 - 1200 AU. 
Assuming a semi-major axis distribution of the form f a (a) oc 
a -1 , our brown dwarf companion frequency of 0.5 ± 0.5% 
(M > 30 Mj) for the range 130 - 520 AU translates to a 
value of 0.83 ± 0.83% for the range 120 - 1200 AU, which 
is in good agreement with the frequency found by McCarthy & 
Zuckerman (2004). This "extrapolated" brown dwarf compan- 
ion frequency may underestimate the true value, if the brown 
dwarf desert does not exist at larger separations (which may 
be the case for low-mass stars in the solar neighbourhood; e.g., 
Gizis et al. 2001). 

6.3. The origin of the brown dwarf desert 

Most stars are formed and reside in binary or multiple stel- 
lar systems. Knowledge about binary and multiple systems in 
young stellar groupings is of fundamental importance for our 
understanding of the star formation process. The formation 
of brown dwarfs and the dearth of brown dwarf companions 
has attained much interest over the last decade. Theories have 
been developed, explaining the existence of the brown dwarf 
desert using migration (Armitage & Bonnell 2002, most ef- 
fective at a ^ 5 AU) or ejection (Reipurth & Clarke 2001) 
of brown dwarfs. The most popular theory that explains the 
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Fig. 9. The companion mass distribution f Ml (M2) for a simulated association. In each panel we show the distribution f Ml (M2) of an association 
consisting of 50 000 binaries for which the primary mass distribution is given by Equation 7. From left to right, the panels show the distribution 
of companion mass in all binaries with stellar primaries, for those of G and K primaries, and for those of A and B primaries, respectively. In 
the top panels we adopted a mass ratio distribution f q {q) oc q~ 033 with < q < 1. In the bottom panels we adopt the same distribution, but with 
the additional constraint that M2 > 0.02 M Q . The brown dwarf regime is indicated with the dashed lines. This figure shows that, even though 
the mass ratio distribution is strongly peaked to low values of q, the substellar-to-stellar companion ratio among intermediate- and high-mass 
stars is very low. 



brown dwarf desert is the embryo ejection scenario. This sce- 
nario predicts ejection of brown dwarfs soon after their for- 
mation by dynamical interactions in unstable multiple systems 
(Reipurth & Clarke 2001). Hydrodynamical calculations (e.g., 
Bate et al. 2003) suggest that star formation is a highly dynamic 
and chaotic process. Brown dwarfs are ejected during or soon 
after their formation. In this theory brown dwarfs can be seen 
as failed stellar companions. 

In the section above we have shown that the small number 
of brown dwarfs among A and late-B members of Sco OB2 can 
be explained with an extrapolation of the mass ratio distribution 
for stellar companions into the brown dwarf regime. There is 
thus no need for a mechanism to remove brown dwarfs. On the 
other hand, the embryo ejection scenario predicts that (at least 
a fraction of) the free-floating brown dwarfs in Sco OB2 have 
been formed as companions to association members. Below, 
we study the consequences in the case that embryo ejection 
has affected the binary population, making use of the detec- 
tion of 28 free-floating brown dwarfs in Upper Scorpius by 
Martin et al. (2004). Under the assumption that this is what 
happened, we roughly estimate the number of primordial bina- 



ries with an A or late-B primary and a brown dwarf companion. 
For comparison between model predictions and observations 
we consider only those companions with an angular separation 
1" < p < 4" and K s < 14 mag, for which our ADONIS and 
NACO observations are complete. 

Martin et al. (2004) present a sample of 104 candidate very 
low mass members, based on DENIS UK photometry, in a re- 
gion of 60 square degrees in US. The authors report spectro- 
scopic observations of 40 of these candidates and show that 
28 are indeed strong candidate members of the US subgroup. 
Under the assumption that Martin et al. (2004) randomly se- 
lected their 40 observed targets out of the sample of 104 can- 
didates, we estimate the total number of brown dwarfs in US 
to be 73 ± 14 in the 60 square degrees region in US. The pro- 
jected area of the three subgroups of Sco OB2 is approximately 
960 square degrees, which gives us an estimate of 1 165 ± 221 
free-floating brown dwarfs in Sco OB2. 

For our estimate of the number of primordial binaries we 
assume that all free-floating brown dwarfs were once compan- 
ions. The number of systems with an A or late-B primary con- 
stitutes a small fraction F m of the total number of binaries, 
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# 


F v ,. 4 „ (%) 


F m (%) 


-^BD.IM.i 


-^BD,IM,i,Iotal 


1 


10.9 


4.93 


6.26 ± 1.19 


7.93 ± 2.05 


2 


10.9 


5.34 


6.78 ± 1.29 


8.45 ±2.11 


3 


10.9 


6.33 


8.04 ± 1.52 


9.71 ± 2.26 


4 


10.9 


0.40 


0.51 ±0.10 


2.18 ± 1.67 


5 


10.9 


2.45 


3.11 ±0.59 


4.78 ± 1.77 


6 


10.9 


2.70 


3.43 ± 0.65 


5.10 ± 1.79 


7 


10.9 


3.19 


4.05 ± 0.77 


5.72 ± 1.84 


8 


10.9 


0.20 


0.25 ± 0.05 


1.92 ± 1.67 



Table 10. An estimate of the number of primordial binaries in 
Sco OB2 with A and late-B primaries and brown dwarf companions 
(with 0.02 M Q < M < 0.08 M Q ) in the angular separation range 
1" < p < 4". Columns 1 and 2 list the model number (cf. Table 8), 
and the fraction of binaries with angular separation 1" < p < 4" (as- 
suming Opik's law). Column 3 lists the fraction F m of primaries in 
the simulated association that is of type A or late-B. Column 4 shows 
the contribution of brown dwarf companions in Sco OB2 with A and 
late-B primaries in the angular separation range 1" < p < 4", inferred 
from the 28 free-floating brown dwarfs in Sco OB2 found by Martin 
et al. (2004), assuming that all brown dwarfs were formed as com- 
panions. Column 5 shows the total number of primordial brown dwarf 
companions with 1" < p < 4" and A or late-B primaries, with the ob- 
served brown dwarf companions (corrected for unseen brown dwarfs 
with K > 14 mag) included. The values of Af B D,iM,i,totai are upper lim- 
its, as it is likely that not all free-floating brown dwarfs were formed 
as companions. 



depending on the mass distribution (see Table 10). Assuming a 
primordial semi-major axis distribution of the form f a (a) oc a -1 
with 15 R < a < 5 x 10 6 R©, about 11% of the brown dwarf 
companions are in the angular separation range 1" < p < 4". 

For each model in Table 10 we calculate how many of the 
1165 ± 221 free-floating have the properties 1" < p < 4" 
and K s < 14 mag, and obtain the contribution A^ B D,iM,i = 
Fi»_4« x Fm x (1 165 ± 221) of the free-floating brown dwarfs 
found by Martin et al. (2004) to the number of primordial bina- 
ries in Sco OB2 with A and late-B primaries and brown dwarf 
companions in the angular separation range 1" < p < 4". We 
estimate the total number of primordial brown dwarf compan- 
ions A^BD,iM,i,totai with 1" < p < 4" of A and late-B primaries 
by adding the observed number of brown dwarf companions, 
corrected for unseen companions with K s > 14 mag (i.e., 
1.67 ± 1.67). 

We have assumed that all free-floating brown dwarfs were 
once companion stars, and therefore obtained upper limits for 

,M,i,totai- By comparing A^bd.im in Table 9 with A^BD,iM,i,totai 
in Table 10 we can derive which primordial mass and mass 
ratio distributions are consistent with the predictions. As not 
necessarily all free-floating brown dwarfs have their origin in a 
binary system, all models with A^bd.im ~ ^V B D,iM,i,totai are con- 
sistent with the predictions (i.e., the current number of brown 
dwarf companions should be less or equal to the primordial 
number of brown dwarf companions). A comparison shows 
that all models are consistent, except model 4 (random pair- 
ing from the Salpeter mass distribution). Under the hypothesis 
that embryo ejection has affected Sco OB2, the current mass 
ratio distribution is slightly shallower than the primordial mass 
ratio distribution, due to the ejection of brown dwarf compan- 



ions. The above derivation gave an estimate of the number of 
primordial binaries with brown dwarf companions, under the 
assumption that the origin of the free-floating brown dwarfs in 
Sco OB2 can be explained with the embryo ejection scenario. 
With our observations we cannot exclude firmly that this hap- 
pened for several of the free-floating brown dwarfs. 

However, if binary formation would result in a mass ratio 
distribution similar to f q (q) oc q~ 33 , the "brown dwarf desert", 
if defined as a deficit of brown dwarf companions relative to 
stellar companions, would be a natural outcome of star for- 
mation. The embryo ejection scenario is not necessary to ex- 
plain the small observed brown dwarf companion fraction in 
this case. 

7. Binarity and multiplicity in Sco OB2 

In Kouwenhoven et al. (2005) we provided a census on bi- 
narity in Sco OB2, consisting of all available data on visual, 
spectroscopic, eclipsing, and astrometric binaries and multi- 
ples. In Table 1 we present an update on the binary statistics in 
Sco OB2. The statistics have been updated with the new results 
presented in this paper, as well as with the binaries recently dis- 
covered by Nitschelm (2004), Jilinski et al. (2006), and Chen 
et al. (2006). 

The multiple system fraction F M , the non-single star frac- 
tion F N s, and companion star fraction Fq are defined as: 

F M = (B + T + ...)/ (S +B + T + ...); (9) 
F NS = (2B + 3T + ...)/ (S + 2B + 3T + ...); (10) 
F c = (B + 2T + ...)/ (S +B + T + ...), (11) 

where S , B, and T denote the number of single systems, binary 
systems and triple systems in the association. In the Sco OB2 
association at least 40% of the systems are multiple. Of the in- 
dividual stars at least 60% is part of a multiple system. Each 
system contains on average Fq ~ 0.5 known companion stars. 
The updated values of F M , F N s, and Fq are slightly larger than 
the values mentioned in Kouwenhoven et al. (2005), respec- 
tively. Note that these frequencies are lower limits due to the 
presence of undiscovered companion stars. 

8. Conclusions 

We have carried out near-infrared JHK S observations of 22 
A and late-B stars in the Sco OB2 association. The obser- 
vations were performed with the NAOS/CONICA system at 
the ESO Very Large Telescope at Paranal, Chile. The observa- 
tions resulted from a follow-up program of our previous work 
(Kouwenhoven et al. 2005), in which we surveyed 199 A and 
late-B Sco OB2 members for binarity with ADONIS. The data 
were obtained with the goal of (1) determining the validity 
of the Ks = 12 mag criterion we used in our ADONIS sur- 
vey to separate companions and background stars, (2) studying 
the behaviour of the companion mass distribution in the low- 
mass regime, and (3) searching for additional companion stars. 
We have included in our analysis the multi-color observations 
of 9 targets observed with ADONIS. In our ADONIS survey, 
these targets were analyzed using their K s magnitude only. The 
main results of our study are: 
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- We detect 72 secondaries around the 3 1 target stars in our 
analysis. By comparing the near-infrared colors with the 
isochrones in the color-magnitude diagram, we find 25 con- 
firmed companion stars, 12 candidate companion stars, and 
35 background stars. 

- For most objects in our ADONIS survey (Kouwenhoven 
et al. 2005) only the K s magnitude was available. We used 
a magnitude criterion to separate companion stars (Ks < 
12 mag) and background stars (K s > 12 mag). With our 
analysis of the 22 NACO targets and 9 ADONIS targets 
with multi-color observations, we estimate the accuracy of 
the Ks = 12 mag criterion. We find that the K s = 12 mag 
criterion is a very useful tool, correctly classifying the sec- 
ondaries in ~ 80% of the cases. 

- We report two candidate brown dwarf companions of 
HIP81972. From their near-infrared magnitudes we infer 
masses of 32 Mj and 63 Mj. The objects are located at an 
angular separation of 7.92" (1500 AU) and 2.79" (520 AU) 
from HIP81972, respectively. Follow-up spectroscopy is 
necessary to determine the true nature of these secondaries. 
Although we are sensitive (but incomplete) to massive plan- 
ets, we classify the faintest secondaries as background stars 
(irrespective of their location in the color-magnitude di- 
agram), because of isochronal uncertainty and the large 
number of faint background stars. 

- In our combined survey of 199 A and late-B members of 
Sco OB2 we detect one confirmed companion star with 
12 mag < K s < 14 mag in the angular separation range 
1" - 4". In this region we detect no other secondary, while 
both the ADONIS and NACO observations are complete. 
This indicates a very low frequency of brown dwarf com- 
panions in the separation range 130-520 AU for late-B and 
A type stars in Sco OB2. 

- Our results are in good agreement with a mass ratio distri- 
bution of the form f q (q) oc q~ 033 . We find a brown dwarf 
companion fraction (for M ^ 30 Mj) of 0.5 ± 0.5% for A 
and late-B stars in Sco OB2. After correction for unseen 
faint companions (M % 30 Mj), we estimate a substellar- 
to-stellar companion ratio of R = 0.06 ± 0.02. 

The number of brown dwarfs among A and late-B members 
of Sco OB2 is consistent with an extrapolation of the (stel- 
lar) companion mass distribution into the brown dwarf regime, 
suggesting that the formation mechanism for stars and brown 
dwarfs is the same. The embryo ejection mechanism does not 
need to be invoked to explain the small number of brown dwarf 
companions among intermediate mass stars in Sco OB2. 
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Appendix A: Results of the NAOS/CONICA survey 

Table A. 1 : Results from our multi-color binarity study among 22 Sco OB2 member stars observed with NACO {top part of 
the table) and the subset of 9 members with multi-color observations in in the ADONIS survey (bottom part of the table). The 
columns show the Hipparcos number (for the targets) and the secondary designation, the /, H, and K s magnitudes, the angular 
separation, and the position angle (measured from North to East). Lower limits to the magnitudes are given if an object is not 
detected in the NACO survey, unless the ADONIS measurement was available (marked with a ★). Entries marked with have 
no available measurement, e.g., because the object is not in the field of view for that filter. For each primary and companion 
star we list the absolute JHKs magnitudes and the mass in columns 7-10. We additionally provide absolute magnitudes and a 
mass estimate for the candidate companions, under the assumption that these are indeed companions. We stress that a significant 
number of the candidate companions may actually be background stars. The 11th column lists the status of the object (p = 
primary, c = confirmed companion star, nc = new confirmed companion star, ? = candidate companion star; b = background 
star). The last column provides additional remarks. A remark "J", "H", or "K" means that the secondary flux in this filter was 
obtained from the image obtained without the NDF, using the PSF from the corresponding image that was obtained with NDF 
(see § 2.4). If the secondary status was obtained without color information, an exclamation mark is placed in the last column. 
The results for the 9 targets with multi-color information in the ADONIS survey are marked with "ADO". 
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> 16.81 


> 16.84 


17.34 


8.13 


36.64 










b 


! 


HIP81972 


5.82 


5.89 


5.87 






-0.56 


-0.49 


-0.51 


4.92 


P 




HIP8 1972-1 


11.63 


10.87 


10.48 


2.02 


313.69 


(5.25) 


(4.49) 


(4.10) 


(0.67) 


? 




HIP81972 -2 


11.30 


10.97 


10.61 


7.02 


258.81 


(4.92) 


(4.59) 


(4.23) 


(0.68) 


? 




HIP81972 -3 


12.54 


11.86 


11.77 


5.04 


213.45 


6.16 


5.48 


5.39 


0.35 


c 


J 


HIP81972 -4 


15.10 


14.43 


13.98 


2.79 


106.94 


8.72 


8.05 


7.60 


0.06 


nc 


JHK 


HIP81972 -5 


16.11 


15.63 


15.26 


7.92 


229.27 


9.73 


9.25 


8.88 


-0.03 


nc 


JHK 


HIP81972 -6 


> 16.58 


16.25 


> 16.61 


8.79 


167.71 










b 


H! 


HIP81972 -7 


> 16.58 


17.12 


> 16.61 


3.58 


33.65 










b 


H! 


HIP81972-8 


> 16.58 


17.28 


> 16.61 


7.44 


265.65 










b 


H! 


HIP83542 


5.34 


4.91 


5.38 






-1.26 


-1.69 


-1.22 


1.10 


P 




HIP83542-1 


** 


9.72 


9.90 


8.86 


196.21 


(**) 


(3.12) 


(3.30) 


(0.91) 


? 




HIP83542 -2 


> 15.54 


15.65 


> 12.13 


9.84 


156.45 










b 


H! 


ADONIS targets with multi-color observations 


HIP53701 


6.30 


6.37 


6.48 






0.79 


0.86 


0.97 


2.84 


P 


ADO 


HIP53701 -1 


9.05 


8.76 


8.86 


3.88 


75.81 










b 


ADO 


HIP53701 -2 


13.06 


12.93 


13.04 


6.57 


120.05 










b 


ADO 


HIP76071 


7.05 


7.10 


7.06 






0.89 


0.94 


0.90 


2.70 


P 


ADO 


HIP76071 -1 


> 11.25 


11.28 


10.87 


0.69 


40.85 


> 5.09 


5.12 


4.71 


0.23 


c 


ADO 


Continued on next page 
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Star 


J 


H 


K s 


P 


PA 


Mj 


M H 


M Ks 


Mass 


Status 


Remarks 




mag 


mag 


mag 


arcsec 


deg 


mag 


mag 


mag 


M Q 






HIP77911 




6.67 


6.71 


6.68 






0.81 


0.85 


0.82 


2.80 


P 


ADO 


HIP77911 


-1 


12.68 


12.20 


11.84 


7.96 


279.25 


6.82 


6.34 


5.98 


0.09 


c 


ADO 


HIP78530 




6.87 


6.92 


6.87 






1.08 


1.13 


1.08 


2.48 


P 


ADO 


HIP78530 


-1 


> 14.50 


14.56 


14.22 


4.54 


139.69 


(> 8.71) 


(8.77) 


(8.43) 


(-0.02) 


? 


ADO 


HIP78809 




7.41 


7.50 


7.51 






1.65 


1.74 


1.75 


2.03 


P 


ADO 


HIP78809 


-1 


11.08 


10.45 


10.26 


1.18 


25.67 


5.32 


4.69 


4.50 


0.30 


c 


ADO 


HIP78956 




7.52 


7.54 


7.57 






1.15 


1.17 


1.20 


2.40 


P 


ADO 


HIP78956 


-1 


9.76 


9.12 


9.04 


1.02 


48.67 


3.39 


2.75 


2.67 


1.16 


c 


ADO 


HIP79124 




7.16 


7.14 


7.13 






1.11 


1.09 


1.08 


2.48 


P 


ADO 


HIP79124 


-1 


11.38 


10.55 


10.38 


1.02 


96.18 


5.33 


4.50 


4.33 


0.33 


c 


ADO 


HIP79156 




7.56 


7.56 


7.61 






1.44 


1.44 


1.49 


2.09 


P 


ADO 


HIP79156 


-1 


11.62 


10.89 


10.77 


0.89 


58.88 


5.50 


4.77 


4.65 


0.27 


c 


ADO 


HIP80238 




7.45 


7.45 


7.34 






1.83 


1.83 


1.72 


1.94 


P 


ADO 


HIP80238 


-1 


7.96 


7.66 


7.49 


1.03 


318.46 


2.34 


2.04 


1.87 


1.67 


c 


ADO 



Table A. 2: Criteria used to determine whether a secondary is a companion star or a background star. Results are listed for secon- 
daries found around the 22 targets observed with NACO (top part of the table) and the 9 targets with multi-color observations in 
the ADONIS dataset (bottom part of the table). Columns 1 and 2 show the secondary designation and the status of the component 
as determined in this paper (c = companion star; ? = candidate companion star; b = background star). Columns 3-5 show the 
compatibility of the location of the object in the color-magnitude diagrams with the isochrones in terms of^ 2 . Confirmed com- 
panions have^ 2 < 2.30 and (confirmed) background stars have^f 2 > 1 1.8. The other secondaries have 2.30 < x 2 < 1 1-8 and are 
labeled "candidate companion". A substantial fraction of these candidate companions may in fact be background stars. Several 
faint (K s > 14 mag) secondaries are only detected in one filter (thus have nox 2 ), and are all assumed to be background stars. 



Star 


Status 


X 1 ]-K s ,Mk s 


Xh-KsMk, 


X]-H,M] 


HIP59502-1 


c 


2.11 


1.09 


0.19 


HIP59502 -2 


b 




75.23 




HIP59502 -3 


b 








HIP60851 -1 


b 


7.77 


1.29 


15.51 


HIP60851 -2 


b 


> 52.81 


0.19 


> 45.54 


HIP60851 -3 


? 




7.54 




HIP60851 -4 


b 




22.34 




HIP60851 -5 


b 




26.37 




HIP60851 -6 


b 








HIP60851 -7 


b 








HIP60851 -8 


b 




237.30 




HIP61265 -1 


? 


4.46 


0.02 


3.83 


HIP61265 -2 


b 


13.77 


1.81 


5.29 


HIP61265 -3 


b 




13.70 




HIP61265 -4 


b 




66.72 




HIP61265 -5 


b 








HIP62026-1 


c 


0.91 


0.09 


0.24 


HIP63204-1 


b 


71.62 


9.64 


17.02 


HIP63204 -2 


c 


1.76 


0.01 


1.38 


HIP67260-1 


c 


0.01 


0.00 


0.01 


HIP67260 -2 


? 




7.20 




HIP67260 -3 


? 


1.45 


3.59 


5.27 


HIP67919-1 


c 


0.02 


0.86 


0.39 


HIP68532-1 


c 


0.93 


1.10 


0.00 


HIP68532 -2 


c 


1.06 


0.02 


1.70 


HIP69113-1 


c 


1.26 


0.09 


0.63 


HIP69113-2 


c 


0.56 


0.04 


1.38 


Continued on next page 
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Star 


Status 


Xj-K s ,M Ks 

ij 


Xh-k s ,m Ks 

j 


Xj-H,M] 


HIP73937-1 


c 




0.01 





HIP73937 -2 


b 


— 


22.37 





HIP78968 -1 


? 


3.88 


1.26 


0.73 


HIP79098 -1 


b 


8.45 


14.45 


39.65 


HIP79410-1 


b 


21.58 


23.51 


19.09 


HIP79739-1 


c 


0.44 


0.09 


0.91 


HIP79771 -1 


c 


1.46 


0.24 


0.52 


HIP79771 -2 


nc 


0.00 


0.04 


0.02 


HIP80142-1 


b 


164.54 


59.81 


33.27 


HIP80142 -2 




— 


— 


2.61 


HIP80474-1 


b 


5.08 


73.31 


36.05 


HIP80799-1 


c 


1.11 


0.12 


0.62 


HIP80896-1 


c 


0.73 


0.02 


0.53 


HIP8 1949-1 


b 


66.30 


7.46 


27.51 


HIP81949 -2 


b 


15.93 


0.91 


8.99 


HIP81949 -3 


b 


> 13.09 


0.02 


> 19.08 


HIP81949 -4 


? 


> 3.88 


5.69 


>5.59 


HIP81949 -5 


b 











HIP81949 -6 


? 


> 9.95 


0.98 


>6.88 


HIP81949 -7 


b 


> 5.16 


> 11.84 





HIP81949 -8 


b 


— 


— 


— 


HIP81949 -9 


b 











HIP81949-10 


b 











HIP81949-11 


b 











HIP81949-12 


b 











HIP8 1972-1 


? 


4.36 


2.27 


0.57 


HIP81972 -2 


? 


1.42 


1.08 


6.03 


HIP81972 -3 


c 


0.56 


1.65 


0.20 


HIP81972 -4 


nc 


0.51 


0.11 


0.13 


HIP81972 -5 


nc 


0.90 


0.07 


0.40 


HIP81972 -6 


b 




> 16.02 




HIP81972 -7 


b 








HIP81972 -8 


b 








HIP83542-1 


? 




8.25 




HIP83542 -2 


b 








ADONIS targets with multi-color observations 


H1P53701 -1 


b 


12.01 


3.72 


0.80 


HIP53701 -2 


b 


33.11 


8.54 


8.16 


HIP76071 -1 


c 




0.57 




HIP77911 -1 


c 


0.73 


0.00 


0.74 


HIP78530-1 


? 




3.10 




HIP78809-1 


c 


0.82 


0.74 


0.02 


HIP78956-1 


c 


0.01 


0.25 


0.33 


HIP79 124-1 


c 


0.22 


0.83 


1.88 


HIP79 156-1 


c 


0.41 


1.94 


0.45 


HIP80238 -1 


c 


0.43 


0.89 


0.06 
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Table A.3: All companion stars identified in our ADONIS and NACO binarity surveys among A and late-B stars in Sco OB2 
(Kouwenhoven et al. 2005, and this paper). The columns show the Hipparcos number of the primary star, the JHKs magnitudes, 
the angular separation, the position angle, the current status of the companion, and the date of observation (dd/mm/yy). If 
measurements are performed in both the ADONIS and NACO surveys, the NACO data are provided. The wide companion 
of HIP77315 at p — 37.37" is HIP77317, another member of Sco OB2. These stars are found to be a common proper motion 
pair (Worley & Douglass 1997), and were both observed in our ADONIS survey. The confirmed and candidate companions for 
which the status is determined using their JHK S photometry, are indicated with "confirmed" and "inconclusive", respectively. 
The candidate companions identified by Kouwenhoven et al. (2005), for which the status is determined using the K s = 12 mag 
criterion, and indicated with "candidate" here. Background stars are not listed here. 



Host primary 


J (mag) 


H (mag) 


K s (mag) 


P(") 


PA (°) 


Companion status 


Date 


HIP50520 


6.39 


2.51 


313.32 


candidate 


06/06/01 


HIP52357 






11.45 


10.04 


72.69 


candidate 


06/06/01 


HIP52357 






7.65 


0.53 


73.01 


candidate 


06/06/01 


HIP56993 


11.88 


1.68 


23.07 


candidate 


06/06/01 


HIP58416 


8.66 


0.58 


166.12 


candidate 


06/06/01 


HIP59413 


8.18 


3.18 


99.83 


candidate 


06/06/01 


HIP59502 


12.35 


11.83 


11.64 


2.94 


26.39 


confirmed 


06/04/04 


HIP60084 


10.10 


46 


329 64 


UclllU ILlclLC 


06/06/01 


HIP60851 

1111 \J\J\J A 




13.63 


13.69 


8.16 


231.46 


inconclusive 

111V-W11V-1 LI J A V V— ' 


06/04/04 

\jyjj \J~j \j~ 


HIP61265 

1111 \J1-£*\J^J 


11.98 


11.66 


11.38 


2.51 


67.15 


inconclusive 


06/04/04 

\J~j \j~ 


HIP61 6^Q 


7.06 


1 87 
1.0/ 


1 89 ztO 


CailUlUaLc 


07/06/01 


HIP61 796 






11.79 


9 89 


108 98 


L tlllU 1 UtllC 


07/06/01 


HIP61796 

i in w a / y \j 






11.86 


12.38 


136.77 


candidate 

V^UllUlUUL^ 


07/06/01 


HIP62002 

1 111 W _ WW — 


7.65 


0.38 


69.24 


candidate 


08/06/01 


HIP62026 

i ill \j 


8.08 


7.90 


7.86 


0.23 


6.34 


confirmed 

VVlllll 1 1 1 L- Ll 


06/04/04 

\J\Jj \J~j \J~ 


HIP62179 

1 111 W_ \- l J 


7.57 


0.23 


282.75 


candidate 

VLll 1L1 1 Ll Ll L L- 


08/06/01 


HIP63204 

1111 vj ^ ' ^ w r 


8.79 


8.51 


8.40 


0.15 


236.56 


confirmed 

VXJlllll 1 1 Ll 


06/04/04 

\Jij \J~ 


HIP6451 5 

i ill w r^y i 


6.94 


0.31 


165.69 


candidate 

LU11U1UULL 


08/06/01 


HIP65822 

1111 \J^J\J£~ _ 


11.08 


1.82 


303.87 


candidate 

LLll 1L1 1 Ll Ll L L. 


08/06/01 


HIP67260 

1111 vj i — WW 




14.04 


14.10 


1.23 


355 65 


inconclusive 


28/04/04 

i-f\jj \J~j \J~ 


HIP67260 


15.84 


14.83 


14.67 


2.33 


77.25 


inconclusive 

lllV-L.'llVl UJl VL 


28/04/04 


HIP67260 


8.88 


8.46 


8.36 


0.42 


229.46 


confirmed 


28/04/04 


HIP67919 


9.98 


9.38 


9.10 


0.69 


296.56 


confirmed 


28/04/04 


HIP68080 


7.19 


1.92 


10.20 


candidate 


05/06/01 


HIP68532 

1111 xJxj^J^J^ 


10.52 


9.85 


9.54 


3.05 


288.50 


confirmed 


28/04/04 

i-t\jj \j~ 


HIP68532 


11.38 


10.94 


10.63 


3.18 


291.92 


confirmed 


28/04/04 


HIP68867 


11.61 


2.16 


284.76 


candidate 


08/06/01 


HIP69113 


10.98 


10.43 


10.29 


5.34 


65.15 


confirmed 


30/04/04 


HIP69113 


11.27 


10.45 


10.30 


5.52 


67.17 


confirmed 


30/04/04 


HIP69749 


11.60 


1.50 


0.84 


candidate 


08/06/01 


HIP70998 


10.83 


1.17 


354.60 


candidate 


06/06/01 


HIP71724 


9.70 


8.66 


23.02 


candidate 


08/06/01 


HIP71727 


7.80 


9.14 


244.96 


candidate 


08/06/01 


HIP72940 


8.57 


3.16 


221.58 


candidate 


06/06/01 


HIP72984 


8.50 


4.71 


260.35 


candidate 


06/06/01 


HIP73937 




8.46 


8.37 


0.24 


190.58 


confirmed 


30/04/04 


HIP74066 


8.43 


1.22 


109.62 


candidate 


08/06/01 


HIP74479 


10.83 


4.65 


154.15 


candidate 


08/06/01 


HIP75056 


11.17 


5.19 


34.51 


candidate 


08/06/01 


HIP75151 


8.09 


5.70 


120.87 


candidate 


08/06/01 


HIP75915 


8.15 


5.60 


229.41 


candidate 


05/06/01 


HIP76001 






7.80 


0.25 


3.17 


candidate 


08/06/01 


HIP76001 






8.20 


1.48 


124.82 


candidate 


08/06/01 


HIP76071 




11.28 


10.87 


0.69 


40.85 


confirmed 


02/06/00, 07/06/01 


Continued on next page 
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Host primary 


J (mag) 


H (mag) K s 


(mag) 


P(") 


PA (°) 


Companion status 


Date 


HIP77315 






7.12 


37.37 


137.32 


candidate 


08/06/01 


HIP77315 






7.92 


0.68 


67.01 


candidate 


05/06/01 


HIP77911 


12.68 


12.20 


11.84 


7.96 


279.25 


confirmed 


02/06/00, 07/06/01 


HIP77939 


8.09 


0.52 


119.13 


candidate 


31/05/00 


HIP78530 




14.56 


14.22 


4.54 


139.69 


inconclusive 


02/06/00, 07/06/01 


HIP78756 


9.52 


8.63 


216.40 


candidate 


02/06/00 


HIP78809 


11.08 


10.45 


10.26 


1.18 


25.67 


confirmed 


03/06/00, 07/06/01 


HIP78847 


11.30 


8.95 


164.02 


candidate 


03/06/00 


HIP78853 


8.45 


1.99 


270.39 


candidate 


08/06/01 


HIP78956 


9.76 


9.12 


9.04 


1.02 


48.67 


confirmed 


03/06/00, 07/06/01 


HIP78968 


14.96 


14.51 


14.26 


2.78 


322.13 


inconclusive 


04/05/04 


HIP79124 


11.38 


10.55 


10.38 


1.02 


96.18 


confirmed 


03/06/00, 07/06/01 


HIP79156 


11.62 


10.89 


10.77 


0.89 


58.88 


confirmed 


03/06/00, 07/06/01 


HIP79250 


10.71 


0.62 


180.92 


candidate 


03/06/00 


HIP79530 


8.34 


1.69 


219.66 


candidate 


31/05/00 


HIP79631 


7.61 


2.94 


127.85 


candidate 


05/06/01 


HIP79739 


12.28 


11.52 


11.23 


0.96 


118.33 


confirmed 


19/06/04 


HIP79771 


12.00 


11.28 


10.89 


3.67 


313.38 


confirmed 


19/06/04 


HIP79771 


12.39 


11.79 


11.42 


0.44 


128.59 


confirmed 


19/06/04 


HIP80142 


16.64 


15.88 




5.88 


119.94 


inconclusive 


04/05/04 


HIP80238 


7.96 


7.66 


7.49 


1.03 


318.46 


confirmed 


02/06/00, 07/06/01 


HIP80324 


7.52 


6.23 


152.46 


candidate 


31/05/00, 03/06/00 


HIP80371 


8.92 


2.73 


140.65 


candidate 


02/06/00, 03/06/00 


HIP80425 


8.63 


0.60 


155.77 


candidate 


08/06/01 


HIP80461 


7.09 


0.27 


285.64 


candidate 


31/05/00 


HIP80799 


10.60 


10.04 


9.80 


2.94 


205.02 


confirmed 


05/05/04 


HIP80896 


11.16 


10.63 


10.33 


2.28 


177.23 


confirmed 


08/06/04 


HIP81624 


7.95 


1.13 


224.28 


candidate 


05/06/01 


HIP81949 




15.62 


14.82 


6.26 


239.37 


inconclusive 


04/05/04, 05/05/04, 08/06/04, 25/06/04 


HIP81949 




15.67 


15.52 


5.27 


340.72 


inconclusive 


04/05/04, 05/05/04, 08/06/04, 25/06/04 


HIP81972 


11.30 


10.97 


10.61 


7.02 


258.81 


inconclusive 


27/06/04 


HIP81972 


11.63 


10.87 


10.48 


2.02 


313.69 


inconclusive 


27/06/04 


HIP81972 


12.54 


11.86 


11.77 


5.04 


213.45 


confirmed 


27/06/04 


HIP81972 


15.10 


14.43 


13.98 


2.79 


106.94 


confirmed 


27/06/04 


HIP81972 


16.11 


15.63 


15.26 


7.92 


229.27 


confirmed 


27/06/04 


HIP83542 




9.72 


9.90 


8.86 


196.21 


inconclusive 


10/09/04 


HIP83693 


9.26 


5.82 


78.35 


candidate 


06/06/01 



